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The  ever  increasing  energy  demands  coupled  with  the  limited  availability  of  fossil  fuels  and  the  detri¬ 
mental  environmental  effects  resulting  from  their  use,  has  guided  research  toward  seeking  alternative 
fuels  to  gradually  substitute  conventional  ones.  Among  these,  biofuels  have  received  increasing  attention 
due  to  their  attractive  features  of  being  renewable  in  nature  and  reducing  the  net  C02  emissions.  Biofuels 
have  been  used  in  conventional  diesel  and  gasoline  engines  either  as  neat  fuels  or  as  supplements. 

Fortunately,  a  relatively  new  combustion  concept  for  internal  combustion  engines,  namely  homo¬ 
geneous  charge  compression  ignition  (HCCI)  combustion,  has  been  evolved  in  parallel  to  the  biofuel 
research.  HCCI  combustion  seems  to  be  able  to  take  advantage  of  the  diverse  properties  of  biofuels,  since 
in  this  combustion  mode  ignition  is  not  externally  instigated,  but  relies  on  the  compression  and  sub¬ 
sequent  autoignition  of  a  fuel-air  mixture.  This  fact  allows  the  utilization  of  different  fuels  or  blends 
thereof,  in  order  to  regulate  the  ignition  point  and  provide  adequate  operation  under  diverse  operating 
conditions. 

This  study  provides  an  overview  of  existing  simulation  models  for  the  simulation  of  biofueled  HCCI 
combustion.  Simulation  models  aid  and  supplement  the  experimental  research  conducted  on  HCCI  com¬ 
bustion,  providing  a  fundamental  insight  into  the  physicochemical  parameters  affecting  performance  and 
emissions  formation.  The  simulation  models  include  single-zone  models,  multi-zone  models,  probability 
based  models,  and  multi-dimensional  models  in  order  of  complexity.  The  vast  majority  of  these  models 
implement  chemical  kinetics  to  simulate  the  combustion  process,  not  only  due  to  the  inherent  depen¬ 
dence  of  HCCI  combustion  on  the  physicochemical  properties  of  the  fuel,  but  also  due  to  the  sometimes 
complex  chemical  structure  of  the  biofuels,  which  include  esters,  ethers  and  alcohols.  The  reaction  paths 
for  these  homologous  series  are  quite  different  from  the  conventional  hydrocarbons  used  to  simulate 
conventional  fuels,  and  provide  the  ground  for  current  and  future  research  work. 

©  2011  Elsevier  Ltd.  All  rights  reserved. 
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Nomenclature 

B 

cylinder  bore  (m) 

h 

molar  specific  enthalpy  (J/kmol) 

i 

any  zone 

k 

species  index 

I< 

number  of  species 

n 

moles 

P 

pressure  (N/m2) 

S 

cylinder  height  (m) 

T 

temperature  (K) 

t 

thickness  (m) 

U 

internal  energy  (J) 

V 

volume  (m3) 

z 

number  of  zones 

Greek  symbols 

X 

air-fuel  equivalence  ratio  (lambda) 

0 

fuel-air  equivalence  ratio  (phi) 

Subscripts 

chem 

chemical 

cyl 

cylinder 

HT 

heat  transfer 

i 

zone  index 

in 

inlet  or  initial 

w 

wall 

Superscripts 

in 

inflow 

out 

outflow 

tr 

transferred 

Abbreviations 

a 

after 

b 

before 

BDC 

bottom  dead  center 

CA 

crank  angle 

CAD 

crank  angle  degrees 

CO 

carbon  monoxide 

CFD 

computational  fluid  dynamics 

Cl 

compression  ignition 

CR 

compression  ratio 

DI 

direct  injection 

DTBP 

di-tertiary  butyl  peroxide 

EDC 

eddy  dissipation  concept 

EGR 

exhaust  gas  recirculation 

EtOH 

ethanol 

EU 

European  union 

EVC 

exhaust  valve  closing 

EVO 

exhaust  valve  opening 

FID 

flame  ionization  detector 

HC 

(unburned)  hydrocarbon(s) 

HCCI 

homogeneous  charge  compression  ignition 

HRR 

heat  release  rate  (J/deg) 

HTR 

high  temperature  reaction 

HTHR 

high  temperature  heat  release 

IA 

injection  angle 

imep 

indicated  mean  effective  pressure  (bar) 

ISFC 

indicated  specific  fuel  consumption  (g/kWh) 

IT 

injection  timing 

IVC 

inlet  valve  closing 

IVO 

inlet  valve  opening 

LES 

large  eddy  simulation 

LTR 

low  temperature  reactions 

LTHR 

low  temperature  heat  release 

MB 

methyl  butanoate 

MD 

methyl  decanoate 

MeOH 

methanol 

MON 

motor  octane  number 

MRG 

methanol  reformed  gas 

MSW 

municipal  solid  waste 

MZ 

multi-zone 

NOx 

nitrogen  oxides 

NVO 

negative  valve  overlap 

OHC 

oxygenated  (unburned)  hydrocarbons 

PDF 

probability  density  function 

PI 

port  injection 

PRF 

primary  reference  fuel 

RES 

renewable  energy  sources 

RON 

research  octane  number 

SRM 

stochastic  reactor  model 

SI 

spark  ignition 

SZ 

single-zone 

TC 

turbocharged 

TDC 

top  dead  center 

toe 

tonnes  oil  equivalent 

UHC 

unburned  hydrocarbons 

1.  Introduction 

During  the  last  decades  two  main  concerns  have  arisen  in  the 
energy  generation  and  transportation  sectors,  namely  the  increas¬ 
ing  energy  demand  coupled  with  the  limited  availability  of  oil 
reserves  and  the  detrimental  environmental  effect  produced  by 
their  exploitation  to  deliver  energy  and  motion.  Nuclear  energy 
could  be  used  to  gradually  displace  fossil  fuel  based  energy  pro¬ 
duction,  but  the  late  accident  in  the  Fukushima  nuclear  power 
plant  has  raised  many  concerns  related  to  the  hazards  of  nuclear 
power  plants;  it  seems  that  the  efforts  toward  the  establishment 
and  the  public  acceptance  of  nuclear  energy  have  yielded  some¬ 
what  to  skepticism.  In  an  effort  to  gradually  displace  fossil  fuels, 
gross  inland  energy  consumption1  of  renewable  energy  sources 
(RES)  has  increased  in  the  European  Union  (27  countries)  in  the 
decade  1999-2009  from  over  5%  to  approximately  9%  [1].  How¬ 
ever,  this  increase  in  the  share  of  renewable  energy  consumption 
seems  to  be  mainly  due  to  the  reduction  in  the  solid  fuel  energy 
consumption  share,  since  the  nuclear  energy  consumption  share 
(~13-14%)  and  the  combined  share  of  petroleum  and  natural  gas 
(~61%)  have  remained  fairly  constant  (Fig.  1). 

The  use  of  fossil  fuels  has  increased  the  energy  dependence  of 
many  countries  around  the  world,  since  oil  production  is  extremely 
localized  with  the  Middle  East  possessing  over  55%  of  the  proved 
oil  reserves  world  wide  as  of  2009  [2,3].  This  impacts  all  countries 
and  citizens  around  the  world,  since  the  oil  prices  are  reflected 
on  virtually  every  commodity.  Moreover,  it  creates  and  sustains 
energy  insecurity  among  countries  with  high  energy  dependence. 

Another  detrimental  effect  of  the  use  of  fossil  fuels  is  related  to 
the  environmental  impact  from  their  use  for  energy  production  and 
transportation.  The  fuel  combustion  in  internal  combustion  engines 
produces  -  among  others  -  C02,  a  major  greenhouse  gas  that  is 


1  Gross  inland  energy  consumption  refers  to  the  total  energy  demand  of  a  country 
or  region.  It  covers  energy  consumption  by  the  energy  sector  itself,  distribution 
and  transformation  losses,  final  energy  consumption  by  end-users  and  statistical 
differences.  It  does  not  include  energy  provided  to  maritime  bunkers  [1  ]. 
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Fig.  1.  Gross  inland  energy  consumption  share  by  fuel  in  the  EU  (27  countries).  Gas: 
almost  exclusively  natural  gas.  Renewable  energy  sources:  biomass  and  renewable 
wastes,  hydro,  geothermal,  wind,  solar,  wave,  tide  and  ocean  energy. 

Data  source:  Eurostat  (online  data  codes:  nrg_100a,  nrglOla,  nrg_102a,  nrg_103a, 
nrg_104a,  nrg_1071a)  [1]. 


known  to  contribute  to  the  global  warming  effect.  The  main  issue 
is  not  the  production  of  C02  itself,  since  the  combustion  of  most 
fuels  leads  to  its  production,  but  the  fact  that  in  the  case  of  fossil 
fuels  the  C02  emitted  cannot  be  reconverted  to  the  original  fuel  by 
any  short-term,  natural  -  and  therefore  energy  neutral  -  process. 

In  view  of  the  first  issues  analyzed,  i.e.  the  energy  dependence 
of  most  of  the  countries,  the  instability  of  the  oil  prices  and  the 
increasing  demand  of  fossil  fuel  energy,  alternative  fuels  have  been 
sought  after.  Among  those,  biofuels  have  received  increasing  atten¬ 
tion  due  to  their  additional  attribute  of  being  renewable.  The  term 
“biofuel”  refers  to  any  fuel  that  can  be  produced  from  biomass, 
such  as  starch,  sugars,  lignocellulosic  material,  animal  fat,  oleagi¬ 
nous  plants,  etc.  The  main  advantage  of  a  biofuel  is  that  it  can  be 
reconstructed  by  natural  processes  using  the  C02  emitted  from  its 
combustion.  This  is  part  of  the  photosynthetic  process  of  the  plants, 
which  use  sunlight  and  C02  to  build  their  cells.  Thus,  biofuels  have 
the  potential  to  be  C02  neutral.  This  is  depicted  in  Fig.  2,  in  which 
a  comparison  is  made  between  the  biofuel  energy  system  and  the 
fossil  fuel  energy  system. 

Initially,  feedstock  such  as  barley,  wheat,  corn,  sugar  beat,  sugar¬ 
cane,  etc.  was  -  and  still  is  -  used  for  ethanol  production;  oleaginous 
plants  such  as  soybeans,  castor  beans,  rapeseed,  Jatropha  and  palm 
were  used  for  biodiesel  production.  However,  the  feedstock  utilized 
for  the  production  of  these  “first  generation”  biofuels  originates 
from  edible  crops.  This  led  to  an  unfavorable  effect  related  to  food 
security  in  terms  of  both  land  use  and  food  prices.  With  over  1  /8th 
of  the  population  worldwide  undernourished  as  estimated  in  2010 
[5,6],  using  part  of  the  edible  crops  to  produce  biofuels  did  not  seem 
a  socially  just  choice.  Added  to  the  land  use  was  the  issue  of  rising 
food  prices  due  to  the  increased  demand,  which  was  instigated  by 
the  alternative  use  of  food  commodities  for  fuel  production  [2]. 

The  food  security  issue  raised  concerns  related  to  the  choice  of 
biofuel  production  routes  and  shifted  research  toward  biofuel  pro¬ 
duction  from  non-edible  feedstock,  mainly  lignocellulosic  material 
or  “plant  biomass”.  The  term  “second  generation”  biofuels  is  used 
to  depict  this  shift  in  the  feedstock  used,  but  not  the  final  fuels 
produced.  The  lignocellulosic  material  is  abundant  in  plants  and 
mainly  consists  of  polysaccharides  [7];  biofuels  can  be  obtained 
from  this  feedstock  by  hydrolysis  and  fermentation,  for  ethanol 


production,  and  by  gasification  to  produce  liquid  hydrocarbons, 
biodiesel,  bioDME  and  bio  Synthetic  Gas.  Additionally,  algae  seem  a 
promising  feedstock  for  the  production  of  second  generation  biofu¬ 
els  and  especially  diesel  fuel.  Their  advantages  are  their  high  growth 
rates,  tolerance  to  diverse  environmental  conditions,  high  C02  trap¬ 
ping  and  fixation  ability,  high  oil  content  -  sometimes  exceeding 
80%  of  their  dry  weight  -  and  the  high  microalgae  yield  per  hectare, 
which  is  almost  an  order  of  magnitude  higher  than  the  palm  oil 
yield  [7].  Owing  to  the  fact  that  other  non-fuel  products  can  also  be 
obtained  from  biomass,  an  integrated  approach  is  being  researched 
to  fully  exploit  the  biomass  available  and  produce  all  possible  prod¬ 
ucts  from  biomass,  including  biofuels  and  a  platform  of  chemicals, 
in  biorefinery  systems  [4,7].  The  production  of  second  generation 
biofuels  is  not  yet  cost  effective  since  there  are  a  lot  of  technical 
barriers  to  be  overcome  [7]. 

In  view  of  the  attractive  attributes  of  biofuels,  efforts  have  been 
made  to  introduce  them  into  the  energy  consumption  systems.  In 
Fig.  3  the  final  energy  consumption2  of  RES  in  EU  is  shown  for  the 
year  2009.  It  can  be  seen  that  the  main  sector  in  which  biofuels  - 
especially  liquid  ones  -  are  currently  used,  is  the  transport  sector 
and  particularly  road  transports.  Fig.  3  also  indicates  that  liquid  bio¬ 
fuels  -  including  bioethanol,  biomethanol,  biodiesel,  bioDME  and 
other  -  are  currently  the  gateway  for  the  penetration  of  RES  into 
the  transport  sector.  The  share  of  liquid  biofuels  in  the  transport 
sector  and  especially  road  transports  has  increased  during  the  last 
decade  as  shown  in  Fig.  4. 

These  data  highlight  the  importance  of  biofuels  and  justify  the 
research  effort  toward  their  utilization  in  conventional  (Cl  and  SI) 
engines  as  well  as  in  Homogeneous  Charge  Compression  Ignition 
(HCCI)  engines.  In  Diesel  engines  fuels  such  as  n-butanol,  ethanol, 
biodiesel,  DME,  either  neat  or  in  blends,  have  been  tested  and/or 
simulated  [8-16].  Some  of  these  biofuels,  such  as  esters  and  DME, 
have  an  inherent  tendency  to  reduce  soot  formation  due  to  their 
chemical  structure;  the  tendency  to  form  soot  precursors  decreases 
as  the  oxygen  content  in  the  fuel  molecule  increases  and/or  the 
number  of  C-C  bonds  decreases  [16].  This  is  due  to  the  fact  that 
the  oxygen  bound  carbon  atoms  tend  to  form  CO  or  C02  rather 
than  participate  in  soot  formation  reactions  [17].  In  spark  ignition 
engines  hydrogen,  biogas,  ethanol,  DME  have  been  tested  and/or 
simulated  [18-21]. 

Apart  from  the  conventional  combustion  modes  a  relatively  new 
concept  has  been  developing  during  the  last  two  decades,  i.e.  homo¬ 
geneous  charge  compression  ignition  combustion  (HCCI).  In  HCCI 
combustion  a  mixture  of  fuel  and  air  is  introduced  into  the  combus¬ 
tion  chamber  early  during  compression  or  during  the  intake  stroke. 
The  mixture  is  compressed  and  autoignites  at  the  point  where  the 
local  thermo-  and  fluid-dynamic  conditions  are  favorable  [22].  The 
main  advantage  of  the  HCCI  concept  relative  to  conventional  com¬ 
bustion  modes  is  the  significant  reduction  of  NOx  emissions  and 
the  absence  of  soot  in  the  exhaust.  Moreover,  the  fast  combustion 
usually  encountered  and  the  unthrottled  operation  are  prone  to 
yield  high  indicated  thermal  efficiency.  However,  HC  and  CO  emis¬ 
sions  are  usually  high  due  to  bulk  quenching  at  low  loads,  and  wall 
quenching  and  crevice  flows  at  higher  loads  [23]. 

Added  to  these  unfavorable  attributes  of  HCCI  is  the  inability 
to  directly  control  the  combustion  phasing  in  HCCI  engines.  This  is 
a  major  obstacle  preventing  the  implementation  of  the  HCCI  con¬ 
cept.  HCCI  ignition  relies  heavily  on  chemical  kinetics  and  largely 
depends  on  the  fuel  used.  In  an  effort  to  control  ignition  and  the 
combustion  rate,  thereby  expanding  the  HCCI  operation  over  a  wide 
range  of  load-engine  speed,  various  means  have  been  used,  such 
as  variable  compression  ratio,  variable  valve  timing,  variable  inlet 


2  Final  energy  consumption  is  the  total  energy  consumed  by  the  end-users,  such 
as  industry,  transports,  households,  services,  etc.  [1  ]. 
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Fig.  2.  Comparison  of  bioenergy  and  fossil  fuel  energy  systems  [4]. 


temperature,  internal  or  external  EGR,  direct  fuel  injection,  etc.  [23  ]. 
A  promising  method  for  the  combustion  control  seems  to  be  the 
use  of  alternative  fuels  and  fuel  blends,  depending  on  the  charac¬ 
teristics  of  the  engine  (compression  ratio,  etc.)  and  the  operating 


conditions  [23-25].  The  main  concept  is  that  fuels  with  different 
autoignition  tendencies  can  be  blended  at  varying  proportions  to 
regulate  the  ignition  point  at  various  load-speed  regions  [26].  This 
is  a  fortunate  evolution  since  the  fuels  used  in  the  blends  may 
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Fig.  3.  Final  energy  consumption  of  RES,  biomass  and  renewable  wastes  by  sector 
in  the  EU  (27  countries)  in  2009.  MSW  (ren.):  municipal  solid  waste  (renewable). 

Data  source:  Eurostat  (online  data  codes:  nrg_1071a,  nrg_1073a)  [1]. 


Fig.  4.  Final  energy  consumption  share  of  liquid  biofuels  in  the  transport  sector  and 
the  road  transport  sector  in  the  EU  (27  countries). 

Data  source:  Eurostat  (online  data  codes:  nrg_100a,  nrg_1073a)  [1]. 
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as  well  be  biofuels.  Thus,  hydrogen  [27],  ethanol  [28-32],  ethers 
[33-37]  and  biodiesel  [38]  have  been  used  in  HCCI  engines,  either 
as  neat  fuels  or  in  blends.  High  cetane  number  fuels  (such  as  DME) 
can  be  used  as  ignition  promoters  in  the  blend  and  high  octane 
number  fuels  (e.g.  ethanol,  methanol,  etc.)  can  be  used  as  ignition 
suppressors.  With  the  use  of  closed  or  open  loop  control  systems, 
the  resulting  fuel  blend  can  provide  the  means  to  control  HCCI  igni¬ 
tion  over  a  wider  range  than  the  one  provided  by  the  individual 
neat  fuels  [25].  This  has  already  been  demonstrated  in  HCCI  engines 
using  blends  of  n-heptane  and  isooctane  [26]. 

As  is  the  case  with  conventional  combustion  types  (SI  and  Cl), 
the  HCCI  combustion  mode  can  be  studied  using  simulation  mod¬ 
els  of  various  degrees  of  complexity.  Simulation  models  describe 
the  fundamental  processes  occurring  within  the  combustion  cham¬ 
ber  and  -  if  accurate  enough  -  can  help  toward  understanding 
the  effects  of  individual  operating  parameters  on  HCCI  combus¬ 
tion.  Such  an  analysis  is  not  always  feasible  in  experiments,  since 
it  is  usually  difficult  to  isolate  the  effect  of  a  single  parameter 
-  be  it  fuel  properties,  inlet  conditions,  etc.  -  without  directly 
or  indirectly  affecting  another  parameter.  Simulation  models  can 
provide  insight  into  the  fundamental  processes  occurring  during 
combustion  and/or  during  the  stage  of  gas  exchange  processes, 
depending  on  the  sophistication  level  of  the  model.  Features  such 
as  combustion  chemistry,  heat  and  mass  transfer  can  be  studied 
so  that  the  crucial  factors  affecting  them  can  be  assessed.  Espe¬ 
cially  in  the  case  of  biofuels,  which  have  not  been  numerically 
studied  extensively  until  recently,  there  is  a  lot  to  be  explored 
as  regards  combustion  chemistry,  emissions  formation,  trans¬ 
port  properties,  etc.  Combustion  chemistry  affects  significantly 
the  mixture  autoignition,  combustion  rate  and  emissions  forma¬ 
tion.  Transport  properties  affect  heat  and  mass  transfer  within  the 
combustion  chamber,  which  determine  the  temperature  field  or 
temperature-mass  distribution  within  the  combustion  chamber. 
Moreover,  in  the  case  of  direct  injection  HCCI,  in  which  the  fuel  is 
directly  injected  into  the  combustion  chamber,  the  fuel  transport 
properties  must  be  adequately  modeled  to  simulate  the  forma¬ 
tion,  break-up,  vaporization  and  mixing  of  the  fuel  jet  with  the 
air.  Such  phenomena  determine  the  mixture  stratification  and  the 
possibility  of  fuel  jet  impingement  on  the  combustion  chamber 
walls. 

HCCI  simulation  models  are  divided  into  single-zone,  multi¬ 
zone  and  computational  fluid  dynamics  (CFD)  models,  in  order  of 
increasing  complexity.  In  single-zone  models  the  main  assumption 
is  the  uniformity  of  the  charge  properties  within  the  combus¬ 
tion  chamber.  Multi-zone  models  consider  multiple  regions  within 
the  combustion  chamber,  each  of  which  is  considered  an  open 
thermodynamic  system  of  uniform  properties.  CFD  models  are 
the  most  sophisticated  ones,  providing  detailed  description  of 
the  in-cylinder  gas  motion  and  thermodynamic  state  on  a  local 
level.  However,  this  increased  level  of  sophistication  is  not  with¬ 
out  compensation,  in  terms  of  increased  computational  time  and 
effort. 

In  the  following  sections  single-zone,  multi-zone  and  CFD  mod¬ 
els,  which  have  been  constructed  for  biofueled  HCCI  combustion 
are  presented.  Since  the  vast  majority  of  HCCI  simulation  models 
use  chemical  mechanisms  of  various  levels  of  detail  to  describe 
the  chemical  kinetics  of  the  ignition  and  combustion  processes,  a 
separate  section  is  devoted  to  the  use  of  chemical  kinetics.  The 
HCCI  simulation  models  presented  in  the  following  sections  are 
included  in  this  study  regardless  of  whether  they  produced  specific 
biofueled-oriented  results.  The  selection  criteria  were  considered 
satisfied  by  the  mere  fact  that  biofueled  HCCI  combustion  was 
simulated.  Upon  consideration  of  the  relatively  short  simulation 
history  of  the  HCCI  engine,  simulation  data  for  biofueled  HCCI 
engines  have  their  own  obvious  merits  and  provide  a  reference  for 
future  studies  nonetheless. 


2.  Chemical  kinetics  modeling  -  reaction  pathways 

Almost  all  of  the  simulation  models  used  in  HCCI  combus¬ 
tion  modeling  and  all  simulation  models  presented  herein  include 
chemical  kinetics.  This  is  the  result  of  the  decisive  role  chemical 
kinetics  play  in  HCCI  combustion.  Since  no  external  means  are  uti¬ 
lized  to  instigate  ignition,  the  latter  relies  on  the  reaction  kinetics  of 
the  fuel.  Ignition  commences  where  the  thermodynamic  conditions 
of  pressure,  temperature  and  species  concentration  are  favorable. 
Moreover,  the  specific  characteristics  of  the  reaction  paths  and 
intermediate  species  also  affect  the  formation  of  pollutants  such  as 
NO*,  CO  and  unburned  HC.  These  issues  provide  the  momentum  for 
the  ongoing  research  oriented  toward  building  chemical  reaction 
paths,  which  are  able  to  describe  in  various  levels  of  detail  the  oxi¬ 
dation  characteristics  of  fuel-air  mixtures  under  various  operating 
conditions.  Compared  to  the  diversity  of  the  biofuels  used,  which 
include  hydrogen,  alcohols,  ethers  and  esters,  the  chemical  mech¬ 
anism  research  has  a  short  history.  The  difficulties  encountered 
involve  the  identification  of  the  primary  and  secondary  reaction 
pathways  as  well  as  the  reduction  of  the  resulting  chemical  mech¬ 
anisms.  This  reduction  is  paramount  if  the  chemical  mechanisms 
are  to  be  used  in  simulation  models,  since  the  simulation  time 
increases  non-linearly  with  the  increase  of  species  and  reactions. 
Such  considerations  are  especially  important  when  using  the  chem¬ 
ical  mechanisms  in  CFD  models,  in  which  the  computational  time 
cost  is  already  aggravated  by  the  solution  of  the  Navier-Stokes 
equations. 

As  regards  the  chemical  mechanisms  used  in  modeling  HCCI 
combustion,  care  should  be  exercised  during  the  choice  of  the  spe¬ 
cific  mechanism,  since  the  various  existing  mechanisms  do  not 
necessarily  provide  similar  results  even  for  well-known  fuels  such 
as  isooctane  [39].  Therefore,  comparisons  of  ignition  timings  and 
global  combustion  rates  must  be  interpreted  with  caution,  espe¬ 
cially  when  differences  between  fuels  are  assessed  [40].  This  also 
holds  true  for  the  comparison  of  NOx  emissions,  which  are  known  to 
be  sensitive  to  peak  temperature  variations  and  consequently  to  the 
ignition  timing  predictions.  Thus,  it  is  essential  for  the  researcher 
to  investigate  the  conditions  under  which  the  mechanisms  were 
developed  and  the  experiments  used  for  their  validation  (shock 
tube,  flow  reactor,  jet  stirred  reactor,  rapid  compression  machine, 
laminar  flame,  HCCI  engine,  etc.).  Preferably,  the  mechanisms  have 
to  be  validated  at  low  equivalence  ratio  and  high  pressure  since 
these  conditions  are  relevant  to  HCCI  combustion  in  engines  [39]. 
Furthermore,  it  has  to  be  kept  in  mind  that  the  fuel  oxidation  chem¬ 
ical  mechanisms  usually  refer  to  surrogates  of  actual  fuels,  which 
cannot  describe  in  exhaustive  detail  all  the  actual  fuel  properties.  It 
is  obvious  that  even  the  detailed  mechanisms  have  their  limitations 
and  their  validation  cannot  be  exhaustive,  considering  the  diverse 
operating  conditions  encountered  in  HCCI  engines  (high  EGR  dilu¬ 
tion,  fuel  injection  in  hot  residuals  during  negative  valve  overlap 
HCCI,  etc.). 

Table  1  includes  chemical  reaction  mechanisms  for  various  bio¬ 
fuels  found  in  the  literature.  Some  of  these  mechanisms  have  not 
been  used  in  the  simulation  studies  presented  herein,  but  are  avail¬ 
able  for  future  research.  However,  the  vast  majority  of  the  oxidation 
mechanisms  included  in  Table  1  has  been  implemented  in  sim¬ 
ulation  models,  which  are  described  in  detail  in  the  subsequent 
sections.  The  citations  describing  each  mechanism  have  been  num¬ 
bered  consecutively,  while  the  citations  of  the  papers  in  which 
these  mechanisms  were  used  have  been  numbered  in  order  of 
appearance  in  the  text. 

3.  Single-zone  modeling 

Single-zone  models  are  often  used  in  HCCI  modeling  due  to 
the  short  computational  time  required  for  the  simulation.  In  the 


Table  1 

Available  reaction  mechanisms  for  biofuel  combustion  simulation. 


#  Number  of  Source  Ref.  Used  in  Ref.a  NO*b  Fuel  Simul.  model 

species/number 
of  reactions 


1 -Butanol  (BuOH) 

1. 

181/1703  and 

Dagaut  et  al. 

- 

573/2701 

[41] 

2. 

234/1399 

Black  et  al.  [42] 

- 

3. 

251/1990 

Dagaut  and 
Togbe [43] 

— 

4. 

1046/4398 

Saisirirat  et  al. 

Saisirirat  et  al. 

BuOH/n-heptane 

SZ 

[44] 

[100] 

Diethyl  ether  (DEE) 

5. 

11 2/484  (+DME, 
EtOH) 

Mack  et  al.  [45] 

Mack  et  al.  [45] 

DEE/EtOH/DTBP 

SZ 

6. 

215/1051 

Yasunaga  et  al. 
[46] 

— 

Dimethyl  ether  (DME) 

7. 

23/23 

Yamada  [47] 

- 

8. 

26/28 

Yao  et  al.  [48], 

Huang  et  al. 

DME 

CFD 

Liang  et  al.  [49] 

[138] 

9. 

27/35  (+MeOH) 

Liang  [50] 

Chenetal.  [82] 

DME-MeOH 

CFD 

30/44 

[133] 

(+MeOH  +  NOx) 

Yao  et  al.  [137] 

DME/MeOH 

CFD 

10. 

28/45 

Kim  et  al.  [51  ] 

- 

11. 

78/336 

Curran  et  al. 

Mack  et  al.  [45] 

[52] 

Yaoetal.  [48] 

DME 

SZ 

Mack  et  al. 

DEE/EtOH 

SZ 

[101] 

Yao  et  al.  [103] 

DME/MeOH 

SZ 

Yao  et  al.  [104] 

- 

dme/ch4 

SZ 

Ogawa  et  al. 

- 

DME/Me0H-H20- 

SZ 

[105] 

EtOH-H2-CH4 

Yamada  et  al. 
[108] 

— 

DME/MeOH-Os 

SZ 

12.  79/351 

Fischer  et  al. 

Hoffman  and 

[83,84] 

DME 

SZ 

[53],  Curran 
et  al.  [54], 

Kaiser  et  al. 

Abraham  [97] 

[55] 

Shudo  and 
Yamada  [106], 
Shudo  et  al. 

- 

dme/h2 

SZ 

[107] 

Kim  et  al.  [120] 

[84] 

DME 

CFD 

Park  [139] 

[84] 

DME 

CFD 

Kong [142] 

[84] 

DME/MeOH 

CFD 

13.  97/457 

Chun-lan  et  al. 

Chun-lan  et  al. 

[85] 

DME 

CFD 

(Detailed)  36/73 
(Reduced) 

[56] 

[56[ 

Engine  typec  CR  Pind  (bar)  Tind  (K)  Load  (0)e  Engine  Validationf  Comments 

speed 

(rpm) 


w  n-heptane 


w  gasoline  sur. 


4S 

16 

353 

0.3 

1500 

QLT 

63%  n-heptane 
37%  BuOH 
(vol.) 

4S 

16.25  1.7 

380-455 

0.3 

1800 

QLT 

4S 

17.0 

1.1 

375 

0.01-0.14 

1400 

QNT 

4S 

17.5 

1.02 

375-380 

0.224,  0.226, 
0.605 

1400 

QNT 

DME  PI,  MeOH 
DI. 

4S 

17.0 

1.16 

375.6 

0.238-0.5 

1400 

QNT 

See  DEE 

4S 

17 

1 

348 

0.215,0.095 

1400 

QNT 

4S 

16.25 

1.7 

- 

0.3,  0.4 

1800 

QLT 

4S 

17 

1 

350 

0DME  ~  0.17 
0MeOH  “ 

0.063 

1400 

QNT 

4S 

17 

1 

348 

0.297 

1400 

No 

4S 

18.8 

1 

293 

0.16,  0.25 

1200 

No 

4S 

8 

1 

378-388 

0.12-0.42 

(0DMe) 

600 

QNT 

4S 

16 

1.07-1.49 

300 

0.4 

2300 

No 

4S 

9.7 

Ambient 

0.27-0.35 

1000 

QLT 

4S 

17.8 

320 

8  mg/stroke 

1500 

QNT 

IT  80°-0°  bTDC 

4S 

17.8 

1-3 

300-460 

0.2-2.0 

1500 

No 

4S 

17.7 

1.1 

340 

0.175-0.435 

960 

QNT 

4S 

— 

1 

350 

0.214 

1500 

QNT 

Emissions 

validation 
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Table  1  ( Continued ) 


# 

Number  of 
species/number 
of  reactions 

Source  Ref. 

Used  in  Ref.a 

NOxb 

Fuel 

Simu 

Ethanol  (EtOH) 

14. 

38/228 

Rohl  and  Peters 

Komninos  and 

[86] 

EtOH 

MZ 

[57] 

Rakopoulos 

[40] 

Komninos  and 

Rakopoulos 

[113] 

[86] 

EtOH 

MZ 

15. 

39/238 

Li  et  al.  [58] 

Sjoberg  and 

Dec  [118] 

— 

— 

— 

16. 

43/235 

Saxena  and 

Viggiano  and 

- 

EtOH 

CFD 

Williams  [59] 

Magi  [145] 

17. 

57/383 

Marinov  [60] 

Szybist  [95] 

Mack  et  al.  [45] 

— 

EtOH 

SZ 

Ng  and 

Thomson  [92] 

[84] 

EtOH 

SZ 

Martinez-Frias 

[84] 

EtOH  (wet) 

SZ 

et  al.  [93] 

Yang  et  al.  [96] 

EtOH 

SZ 

Mack  et  al. 

- 

EtOH/DEE 

SZ 

[101] 

Ogawa  et  al. 
[105] 

Komninos  and 

[86] 

EtOH 

MZ 

Rakopoulos 

[113] 

Sjoberg  and 

Dec  [118] 

Yu  et  al.  [143] 

EtOH 

CFD 

Joelsson  et  al. 
[144] 

— 

EtOH 

CFD 

18. 

58/310 

Curran  et  al. 

Sjoberg  and 

- 

EtOH 

MZ 

[61] 

Dec  [117] 

Sjoberg  and 

Dec  [118] 

- 

EtOH 

MZ 

19. 

112/484 

Curran  in 

Mosbach  et  al. 

- 

EtOH/DEE 

SRM 

Mosbach  et  al. 
[62] 

[62] 

20. 

167/1591, 

Dagaut  and 

- 

564/2589 

Togbe [63] 

21. 

235/1866 

Dagaut  and 
Togbe [64] 

— 

22. 

1046/4398 

Saisirirat  et  al. 

Saisirirat  et  al. 

- 

EtOH/n-heptane 

SZ 

[44] 

[100] 

Ethyl  acetate  (Ethyl  Acet.) 

23. 

23/142 

Gasnot  et  al. 

Cotino  and 

- 

Ethyl  Acet. 

CFD 

[651 

Jeanmart  [150] 

Hydrogen 

24. 

10/19 

Mueller  et  al. 

[66] 

— 

25. 

10/19 

Marinov  et  al. 

Goldsborough 

- 

h2 

SZ 

[67] 

and  Van 

Blarigan  [89] 
Komninos  et  al. 

[86] 

h2 

MZ 

[112] 


cn 

CD 

4^ 


Engine  typec 

CR 

P,„d  (bar) 

Tind  (K) 

Load  (0)e 

Engine 

speed 

(rpm) 

Validationf 

Comments 

4S 

21 

1 

-393 

0.1-0.28 

994 

QNT 

4S 

21 

1 

-393 

0.15-0.26 

994 

QNT 

- 

- 

- 

- 

- 

- 

- 

Examined,  not 
used 

4S 

21 

1 

363-393 

0.16-0.26 

1000 

QNT 

4S 

7-29 

1 

325,500-550 

1,  -0.25 

900,  2000 

No 

NVO 

See  DEE 

4S 

15 

1 

380-660 

0.2-0. 5 

1000 

No 

EGR,  neat  EtOH 
or  Reformed 

4S 

16 

3 

420  (ivc) 

1  (4* EtOH  ) 

1800 

No 

TC 

4S 

13 

— 

— 

— 

1400 

No 

NVO 

4S 

16.25 

1.7 

0.3,  0.4 

1800 

QLT 

See  DME 

4S 

21 

1 

-393 

0.15-0.26 

994 

QNT 

- 

- 

- 

- 

- 

- 

- 

Examined,  not 
used 

4S 

17.2 

1 

354,  400 

0.303 

1200 

QNT 

4S 

17.2,17 

0.9-1 

354,388 

0.303 

1200 

QNT 

4S 

12.8 

0.84-1.81 

418-443 

0.25-0.6 

500-2400 

QNT 

4S 

12.8 

1 

427 

0.4 

1200 

QNT 

4S 

16.25 

1.18 

400-420 

0.3,  0.4 

1800 

QNT 

w  n-heptane 

w  gasoline  surr. 

4S 

16 

- 

353 

0.3 

1500 

QLT 

63%  n-heptane 
37%  EtOH  (vol.) 

4S 

13.3 

1.35  boost 

462 

1.5-4  bar 
imepg 

1200 

QNT 

Validation  with 

iso-octane 

2S 

3-38 

1.5 

300 

0.15-0.6 

- 

No 

Free  piston 
engine 

4S 

18-20 

1-2 

390-410 

0.2-0.4 

1000-3000 

No 
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Table  1  ( Continued ) 


#  Number  of 

Source  Ref. 

Used  in  Ref.a 

NO/ 

Fuel 

Simul.  model 

Engine  typec 

CR 

Pind  (bar) 

Tj  (K) 

Load  (0)e 

Engine 

Validationf 

Comments 

species/number 

speed 

of  reactions 

(rpm) 

26.  11/19 

O’Connaire 
et  al.  [68] 

- 

27.  11/23 

Westbrook  and 

Noda  and 

[83] 

h2 

MZ 

4S 

20 

345-410 

0.1-1.0 

1600 

No 

Dryer [69] 

Foster  [111] 

28.  11/25 

Kee  et  al.  [70] 

Fiveland  and 
Assanis  [88] 

— 

h2 

SZ 

4S 

15 

1 .2.1.5 

400-800 

0.15-0.4 

1500 

No 

29.  12/30 

Li  and  Karim 

Liu  and  Karim 

[71,72] 

h2 

CFD 

4S 

16-20 

0.85 

280-320,515 

0.5 

800,  900 

QLT 

[71],  Liu  [72] 

[149] 

Methanol  (MeOH) 

30.  20/84 

Li  et  al.  [73] 

31.  22/89 

Held  and  Dryer 
[74] 

Yao  et  al.  [103] 

See  DME 

32.  26/84 

Westbrook 

Ogawa  et  al. 

See  DME 

et  al.  [75,76] 

[105] 

33.  52/326 

Lindstedt  and 

Yang  et  al.  [96] 

- 

MeOH 

SZ 

4S 

13 

- 

- 

- 

1400 

No 

NVO,  32 

Meyer  [77] 

Species/167 

Reactions 

Methyl  butanoate  (MB) 

34.  41/150  (neat) 

Brakora  et  al. 

Brakora  and 

[87] 

MB  +  n-heptane 

SZ 

4S 

18 

1.1 

300-400 

0.298-1.0 

1200 

No 

56  Species/169 

53/156  (w 
n-heptane) 

[78] 

Reitz  [98] 

Reactions 

Um  and  Park 
[152] 

[84] 

MB  +  n-heptane 

CFD 

4S 

17.8 

1-3  boost 

300-460 

0.2-2.0 

1500 

No 

35.  275/1545 

Dooley  et  al. 

[79] 

— 

36.  295/1498 

Gai'l  et  al.  [80] 

Szybist  et  al. 

[38] 

[84] 

MB  +  n-heptane 

SZ 

4S 

10.5 

1 

433-498 

0.38-0.48 

1800 

QNT 

Tin  variable 

Methyl  decanoate  (MD) 

37.  3036/8555 

Herbinet  et  al. 

Hofmann  and 

[83,84] 

MD 

SZ 

4S 

16 

1.16-1.65 

310,315 

0.4 

1670 

No 

[81] 

Abraham  [97] 

QLT,  qualitative  validation-mostly  validation  of  predicted  trends;  QNT,  quantitative  validation-includes  direct  comparison  of  simulated  quantities  to  experimental  measurements;  NUM,  numerical  validation-comparison  of  the 
HCCI  simulation  results  using  a  reduced  mechanism  against  predictions  from  a  detailed  mechanism. 

a  The  source  mechanisms  may  have  been  modified  for  their  use  in  each  study.  Moreover,  some  of  the  mechanisms  may  have  been  updated  since  the  publication  of  each  study.  Thus,  the  number  of  species  and  reactions  are 
subject  to  change.  Symbol  used  in  this  column  indicates  mechanisms,  which  are  not  used  in  the  studies  presented  herein. 
b  Symbol  in  NO*:  information  not  available  in  the  study. 
c  4S:  four  stroke  engine,  2S:  two  stroke  engine 
d  Subscript  “in”  refers  to  inlet  or  initial  conditions 

e  Load  is  expressed  in  terms  of  the  fuel-air  equivalence  ratio,  unless  explicitly  stated  otherwise. 
f  Validation  refers  to  the  corresponding  HCCI  simulation  model  and  not  to  the  reaction  mechanism(s). 
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Fig.  5.  Schematic  of  the  First  Law  of  Thermodynamics  for  single-zone  models. 


single-zone  (or  zero-dimensional)  model  the  engine  cycle  is  mod¬ 
eled  using  the  First  Law  of  Thermodynamics,  taking  into  account 
the  work  absorbed/produced  by  the  piston  and  -  in  some  cases  - 
the  heat  exchanged  between  the  charge  and  the  combustion  cham¬ 
ber  walls,  as  shown  in  Fig.  5.  Mass  blow-by  is  usually  neglected  as 
shown  in  the  figure. 

The  thermodynamic  (pressure,  temperature,  composition,  etc.) 
and  transport  properties  (viscosity,  conductivity,  etc.)  of  the  mix¬ 
ture  are  considered  uniform.  Since  any  temperature  stratification 
within  the  combustion  chamber  is  neglected,  the  single-zone 
models  tend  to  overestimate  the  global  heat  release  rate  and,  con¬ 
sequently,  the  pressure  rise  rate  and  the  maximum  combustion 
pressure.  Moreover,  since  the  mixture  composition  is  also  consid¬ 
ered  uniform,  the  single-zone  model  is  incapable  of  describing  the 
HC  and  CO  formation  processes,  unless  these  are  due  to  global  par¬ 
tial  combustion  of  the  mixture;  under  normal  HCCI  combustion 
conditions,  the  formation  of  these  pollutants  is  believed  to  occur  in 
the  relatively  low  temperature  regions  of  the  combustion  chamber 
as  a  result  of  mass  transfer  between  cold  and  hot  regions,  in  which 
any  HC  transferred  are  partially  burned  producing  CO  [23]. 

Apart  from  these  inherent  limitations,  single-zone  models  can 
be  used  to  provide  an  estimation  of  the  ignition  timing  and  the 
formation  of  NOx  emissions.  This  is  due  to  the  fact  that  both  of 
these  phenomena  are  primarily  controlled  by  the  temperature  of 
the  hottest  regions  of  the  combustion  chamber,  which  is  close  to 
the  mean  gas  temperature  of  the  mixture.  Moreover,  single-zone 
models  can  be  used  to  assess  the  relative  importance  of  chemical 
reactions  and  dominant  species  that  are  part  of  a  chemical  reaction 
mechanism  describing  the  fuel-air  oxidation.  Thus,  detailed  reac¬ 
tion  mechanisms  can  be  compared  to  reduced  ones  under  the  same 
initial  conditions  and  assumptions  provided  by  the  single-zone 
simulations.  The  latter  features  are  especially  important  in  the  case 
of  biofuels,  the  combustion  chemistry  of  which  is  under  continuous 
investigation.  The  relevant  studies  are  perplexed  by  the  variety  of 
the  biofuels,  which  include  organic  compounds  of  diverse  homolo¬ 
gous  series,  such  as  alcohols,  ethers,  and  esters.  This  is  a  significant 
departure  from  the  combustion  chemistry  of  the  common  pure 
hydrocarbons,  which  are  usually  used  as  surrogates  for  the  simu¬ 
lation  of  diesel  or  gasoline  combustion  and  do  not  include  oxygen 


atoms  in  their  molecules.  Besides  the  study  of  the  combustion 
chemistry  of  neat  biofuels,  single-zone  models  can  be  used  to  study 
the  characteristics  of  various  fuel  blends  or  the  effect  of  additives 
on  biofuel  HCCI  combustion.  This  is  an  important  feature,  since  the 
use  of  fuel  blends  and  additives  may  provide  the  means  for  ignition 
timing  and  combustion  rate  control.  Issues  like  the  aforementioned 
provide  the  impetus  for  the  ongoing  HCCI  research  [88]. 

3.1.  Hydrogen 

Goldsborough  and  Van  Blarigan  [89]  studied  numerically  a  free- 
piston  hydrogen-fueled  HCCI  engine.  The  researchers  modeled 
piston  dynamics  and  the  gas  exchange  process  in  their  free-piston, 
two-stroke  engine.  They  also  used  the  HCT  [90]  chemical  kinetics 
software  to  simulate  a  zero-dimensional  thermodynamic  process 
throughout  the  engine  closed  cycle.  The  hydrogen  oxidation  was 
accounted  for  using  the  chemical  mechanism  of  Marinov  et  al. 
[67].  Heat  transfer  was  accounted  for  by  using  Woschni’s  corre¬ 
lation  [91],  in  which  the  heat  transfer  constants  were  modified 
to  match  experimental  data  from  hydrogen  combustion  in  a  rapid 
compression  machine.  It  was  found  that  the  engine  operation  was 
significantly  dependent  on  the  scavenging  process,  which  altered 
the  amount  of  fuel-air  mixture  trapped  within  the  combustion 
chamber,  the  initial  temperature  of  the  charge,  and  the  amount 
of  fresh  charge  that  was  expelled  directly  through  the  exhaust  port 
during  this  scavenging  process.  Overall,  the  engine  was  found  to 
produce  low  NOx  emissions  and  to  provide  almost  constant  volume 
combustion.  The  rapid  combustion,  which  resulted  in  the  simula¬ 
tions,  could  be  partly  attributed  to  the  nature  of  the  single-zone 
model  used. 

Fiveland  and  Assanis  [88]  constructed  a  zero-dimensional 
numerical  model  for  the  simulation  of  the  complete  engine  cycle 
using  methane  or  hydrogen  as  a  fuel.  The  hydrogen  combustion 
chemistry  scheme  included  11  species  and  25  chemical  reactions, 
and  was  modeled  using  the  CHEMKIN  libraries  [70]  adapted  for  a 
variable  volume  plenum.  Heat  transfer  was  calculated  using  a  sub¬ 
set  of  the  k-s  model  under  the  isotropic  turbulence  assumption. 
This  method  yielded  the  characteristic  velocity  that  was  used  to 
calculate  the  heat  transfer  coefficient. 

Simulation  results  for  a  hydrogen-fueled  four-stroke  engine 
with  a  compression  ratio  of  1 5  and  1 .5  bar  intake  pressure  required 
425 1<  intake  charge  temperature  to  produce  ignition  nearTDC.  This 
high  intake  temperature  requirement,  at  a  moderate  compression 
ratio  and  at  a  low  fuel-air  equivalence  ratio  of  0.3,  was  necessary 
since  hydrogen  ignition  occurred  at  about  1050-1100  K  in  the  sim¬ 
ulation.  Under  these  conditions  the  researchers  identified  peaks  in 
H,  OH  and  H02  fractions  near  the  ignition  point.  The  simulation 
model  was  also  used  in  a  parametric  investigation,  in  which  it  was 
found  that  decreasing  the  initial  temperature  or  the  fuel-air  equiva¬ 
lence  ratio,  beyond  a  certain  point,  retarded  ignition  and  eventually 
led  to  misfire. 

3.2.  Alcohols 

Ng  and  Thomson  [92]  constructed  a  single-zone  model  to  sim¬ 
ulate  ethanol  HCCI  combustion  and  especially  the  effect  of  ethanol 
reformed  gas  on  the  ignition  requirement  and  NOx  emissions.  They 
incorporated  Marinov’s  set  of  chemical  reactions  [60]  to  simulate 
ethanol  oxidation  and  integrated  a  NOx  reaction  mechanism  from 
GRI-Mech  [84]  to  simulate  NOx  formation.  For  the  estimation  of 
heat  losses  they  used  a  convective  heat  transfer  correlation  (similar 
to  Newton’s  Law  of  cooling). 

The  researchers  conducted  a  parametric  study,  in  which  the 
parameters  varied  were  the  fuel-air  equivalence  ratio,  the  EGR 
amount  and  the  degree  of  reforming,  i.e.  the  degree  of  ethanol  con¬ 
version  into  CO  and  H2.  Thus,  the  initial  charge  contained  ethanol, 
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CO,  H2  and  EGR  at  various  proportions.  It  was  found  that  fuel 
reforming  resulted  in  an  expansion  of  the  HCCI  operating  region, 
the  boundaries  of  which  were  determined  by  misfire  and  unac¬ 
ceptable  NOx  emissions.  However,  reforming  tended  to  increase 
NO*  emissions  due  to  the  resulting  high  peak  combustion  temper¬ 
ature.  In  relation  to  reformed  gas,  EGR  was  more  effective  toward 
expanding  the  simulated  operating  range  at  moderate  equivalence 
ratios  (0.2-0.5). 

Martinez-Frias  et  al.  [93]  investigated  the  direct  use  of  wet 
ethanol  using  a  simulation  model.  The  use  of  wet  ethanol  as  a  fuel 
reduces  the  energy  demands  for  its  production,  since  the  distilla¬ 
tion  and/or  dehydration  processes,  which  are  used  to  remove  the 
water  and  provide  neat  ethanol,  require  a  significant  amount  of 
energy  (37%  of  the  total  energy  in  ethanol  products  for  corn-based 
ethanol  production).  The  objective  of  the  study  was  to  determine 
whether  benefits  could  be  obtained  by  the  direct  utilization  of  wet 
(as  opposed  to  neat)  ethanol  into  a  turbocharged  HCCI  engine.  The 
researchers  constructed  a  simulation  model  for  a  TC  HCCI  engine 
fueled  with  wet  ethanol  (35vol.%  EtOH  in  H20).  Apart  from  the 
TC,  the  full  engine  system  model  also  incorporated  an  exhaust  gas 
regenerator,  a  fuel  vaporizer  and  a  catalytic  converter. 

The  HCCI  engine  was  modeled  via  a  single-zone  model,  using 
CHEMKIN  [94]  to  solve  for  chemical  kinetics.  The  ethanol  reac¬ 
tion  mechanism  was  that  of  Marinov  [60]  and  the  NOx-related 
reactions  were  taken  from  GRI-Mech  [84].  Heat  transfer  was  also 
incorporated  into  the  reciprocating  engine  model  assuming  a  wall 
temperature  of  430  K.  Using  this  model,  the  researchers  found  that 
a  minimum  ethanol-in-water  fraction  of  35  vol.%  (liq.)  was  neces¬ 
sary  for  a  sustainable  simulated  HCCI  operation.  This  provided  an 
overall  brake  efficiency  of  38.7%  for  the  entire  simulation  system. 
The  researchers  estimated  that  this  direct  use  of  35%  wet  ethanol 
significantly  reduced  the  energy  required  for  fuel  production  -  rel¬ 
ative  to  neat  ethanol  production  -  by  eliminating  the  need  for 
dehydration  and  limiting  the  distillation  needs  to  3%  of  the  total 
energy  in  ethanol  products,  in  the  case  of  corn-based  ethanol  pro¬ 
duction. 

In  Ref.  [95]  Szybist  used  CHEMKIN-PRO  [94]  to  simulate  the  HCCI 
combustion  of  n-heptane,  isooctane,  toluene  and  ethanol.  The  sim¬ 
ulation  model  was  switched  to  the  adiabatic  mode.  The  mechanism 
used  for  the  simulation  of  ethanol  combustion  was  that  of  Mari¬ 
nov  [60]  including  56  species  and  351  reactions.  The  author  used 
this  model  to  simulate  HCCI  combustion  with  the  aforementioned 
fuels  under  normal  naturally  aspirated  conditions  with  air  dilution 
and  under  negative  valve  overlap  (NVO)  conditions.  In  the  latter 
case,  the  mixture  had  a  higher  temperature  at  BDC  (>500  K)  and  a 
significant  amount  of  diluents  such  as  02,  C02,  H20  and  N2.  These 
conditions  resemble  HCCI  under  NVO  operation,  in  which  a  portion 
of  the  exhaust  gases  is  trapped  within  the  combustion  chamber  and 
raises  the  temperature  of  the  incoming  fuel-air  charge.  The  results 
of  the  study  showed  that  with  low  intake  temperature  (325  K)  and 
no  EGR,  the  most  resistant  fuel  to  ignite  was  ethanol  requiring  the 
highest  compression  ratio.  However,  under  NVO  conditions  toluene 
was  the  most  difficult  to  ignite.  Moreover,  with  NVO  the  igni¬ 
tion  timing  between  fuels  converged  to  within  a  few  CA  degrees, 
with  the  exception  of  n-heptane  at  high  02  dilution.  On  the  con¬ 
trary,  ethanol  was  essentially  insensitive  to  the  02  dilution,  under 
the  conditions  of  the  simulation,  indicating  a  zero-order  reaction 
rate  with  respect  to  02.  The  researcher  concluded  that  ignition 
in  HCCI  combustion  under  NVO  conditions  presented  distinct  fea¬ 
tures,  which  were  not  entirely  consistent  with  the  RON  and  MON 
parameters. 

Yang  et  al.  [96]  also  investigated  NVO  HCCI  combustion.  They 
used  a  single-zone  adiabatic  model  to  simulate  the  reactions  tak¬ 
ing  place  within  the  combustion  chamber  with  the  injection  of 
gasoline,  ethanol  or  methanol  during  the  NVO  period.  The  reac¬ 
tion  kinetics  was  incorporated  via  the  mechanisms  of  Marinov 


[60]  for  ethanol,  and  Lindstedt  and  Meyer  [77]  for  methanol.  The 
researchers  studied  both  the  thermal  and  the  chemical  effect  of 
the  fuel  injection  during  the  NVO  period.  They  found  that  during 
this  period  both  heat  and  active  species  such  as  H2  and  CO  can 
be  produced,  depending  on  the  fuel  type,  oxygen  availability  and 
the  in-cylinder  temperature  of  the  residuals  during  injection.  At 
relatively  high  oxygen  availability  in  the  residuals  the  tendency 
was  to  produce  more  heat  than  active  species,  since  the  fuel  was 
oxidized.  In  oxygen-poor  residuals  only  part  of  the  fuel  was  oxi¬ 
dized  producing  active  species.  Methanol  and  ethanol  were  less 
sensitive  than  gasoline  to  02  availability  owing  to  their  oxygenated 
nature.  Methanol  showed  the  greatest  tendency  in  forming  active 
species. 

3.3.  Ethers 

In  [48]  Yao  et  al.  used  the  SENKIN  code  of  the  CHEMKIN  pack¬ 
age  to  simulate  DME  combustion.  The  single-zone  model  used 
was  considered  adiabatic  and  the  chemical  mechanism,  which  was 
used  to  describe  the  DME  oxidation,  was  that  of  Lawrence  Liver¬ 
more  National  Laboratory  including  79  species  and  336  reactions 
[52-55].  The  scope  of  the  study  was  to  analyze  the  reaction  path¬ 
ways  of  DME  oxidation.  The  researchers  validated  the  simulation 
results  against  experimental  data,  which  indicated  that  the  com¬ 
bustion  rate  and  the  peak  combustion  pressure  were  overestimated 
due  to  the  nature  of  the  single-zone  model.  However,  the  ignition 
timing  was  adequately  predicted,  implying  a  sufficient  descrip¬ 
tion  of  the  reaction  pathways.  It  was  found  that  DME  presented  a 
two  stage  auto-ignition  with  low  temperature  reactions  (LTR)  fol¬ 
lowed  by  high  temperature  reactions  (HTR).  HCHO  and  OH  were 
the  primary  radicals  during  LTR.  DME  ignition  was  controlled  by 
the  H202  decomposition  to  two  OH  radicals,  which  provided  a 
chain  branching  step.  These  hydroxyls  were  the  primary  attack¬ 
ers  for  the  H  atom  abstraction  from  DME  and  the  oxidation  of 
CO.  At  relatively  leaner  DME  concentration  the  LTR  was  not  sig¬ 
nificantly  affected  but  the  CO  oxidation  deteriorated,  indicating  a 
temperature  sensitivity  of  the  CO  oxidation  process.  This  led  to  the 
generally  accepted  conclusion  that  CO  emissions  in  HCCI  engines 
originate  from  the  relatively  cold  and/or  lean  regions  of  the  com¬ 
bustion  chamber,  such  as  the  crevices  and  the  vicinity  of  the  wall, 
in  which  the  local  temperature  is  adequately  high  to  induce  partial 
oxidation  of  the  existing  HC  but  not  high  enough  to  consume  the 
resulting  CO. 

Hofmann  and  Abraham  [97]  also  used  the  internal  combustion 
engine  model  included  in  the  CHEMKIN  package  [94]  to  simulate 
n-heptane,  DME  and  methyl  decanoate  HCCI  combustion.  The  lat¬ 
ter  was  studied  as  a  biodiesel  surrogate.  The  simulation  model 
incorporated  Woschni’s  correlation  [91]  for  the  estimation  of  heat 
losses  from  the  combustion  chamber.  The  primary  objective  of  the 
study  was  the  investigation  of  ignition  delay  and  NO  formation 
for  the  three  fuels.  For  this  purpose,  the  researchers  used  Her- 
binet  et  al.  [81]  detailed  chemical  kinetic  mechanism  for  methyl 
decanoate  -  a  biodiesel  surrogate  -  consisting  of  3036  species 
and  8555  reactions,  and  Livermore  National  Laboratory  scheme 
for  DME  oxidation  including  79  species  and  351  reactions  [53-55]. 
NO  formation  was  simulated  using  two  alternative  schemes;  the 
extended  reaction  mechanism  of  Zeldovich  [83  ]  and  the  NO  kinetics 
from  GRI-Mech  v.3  mechanism  [84].  The  Zeldovich  NO  formation 
scheme  provided  reduced  NO  concentration  relative  to  the  more 
extensive  NO  scheme  of  GRI-Mech  for  all  fuels.  Moreover,  it  was 
demonstrated  that  DME  showed  the  shortest  ignition  delay  fol¬ 
lowed  by  methyl  decanoate  and  n-heptane.  Thus,  DME  required 
higher  engine  speed  and  the  lowest  initial  temperature.  Although 
the  primary  consideration  for  NO  formation  was  the  peak  combus¬ 
tion  temperature,  fuel  chemistry  was  also  found  to  play  a  secondary 
role. 
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3.4.  Biodiesel 

Apart  from  Hofmann  and  Abraham  mentioned  in  the  preced¬ 
ing  sub-section,  biodiesel  HCCI  combustion  was  also  simulated  by 
Brakora  and  Reitz  [98]  using  the  SENKIN  code  of  the  CHEMKIN 
package  [94].  The  single-zone  model  was  considered  adiabatic  for 
the  purpose  of  the  study,  which  was  to  investigate  the  formation 
of  NOx  from  biodiesel-  and  diesel-fueled  HCCI  engines.  N-heptane 
was  used  as  the  diesel  surrogate,  and  a  mixture  of  2/3  n-heptane 
and  1/3  methyl  butanoate  was  used  to  simulate  biodiesel  combus¬ 
tion.  This  proportion  of  n-heptane  and  methyl  butanoate  in  the 
biodiesel  surrogate  was  found  to  yield  a  better  approximation  for 
the  actual  C:H:0  ratio  found  in  the  real  biodiesel,  relative  to  neat 
methyl  butanoate  or  methyl  decanoate.  The  latter  two,  being  rela¬ 
tively  short  chain  esters,  provided  a  much  higher  C:0  ratio  than 
the  actual  biodiesel  fuel.  The  chemical  mechanism  used  for  the 
biodiesel  surrogate  consisted  of  56  species  and  169  reactions  [78]. 
NO*  formation  was  incorporated  via  a  7  species  and  19  reactions 
scheme,  which  was  developed  by  Yoshikawa  and  Reitz  [87]  relying 
on  reactions  from  the  GRI-Mech  mechanism  [84]. 

Brakora  and  Reitz  found  that  the  ignition  timing  played  a  cru¬ 
cial  role  in  the  resulting  NO*  emissions.  With  an  equivalent  energy 
release  and  for  the  same  ignition  timing,  the  two  surrogates  pro¬ 
vided  roughly  the  same  peak  combustion  temperatures  (within 
approximately  20  K)  and  yielded  similar  NO*  emissions  at  EVO. 
There  was  a  difficulty  in  connecting  the  oxygen  content  of  the  fuel 
to  the  produced  NOx  emissions,  especially  at  lean  mixtures.  Under 
stoichiometric  conditions,  the  biodiesel  surrogate  yielded  a  26% 
increase  in  NOx  relative  to  the  diesel  surrogate.  However,  the  pre¬ 
dicted  difference  in  peak  combustion  temperature  was  about  13  K 
lower  for  biodiesel,  suggesting  that  the  fuel-bound  oxygen  could 
play  a  role  in  the  formation  of  NOx  under  limited  02  concentrations. 
It  was  also  found  that  under  the  same  energy  content,  ignition  tim¬ 
ing  and  total  N:0  ratio  in  the  mixture,  for  the  two  fuels,  the  diesel 
surrogate  yielded  higher  NOx  emissions.  This  suggested  that  hav¬ 
ing  the  same  amount  of  02  available  in  the  air,  as  opposed  to  being 
fuel-bound,  increased  the  NOx  emissions  due  to  the  alteration  of 
the  heat  capacity  of  the  mixture,  which  resulted  in  higher  flame 
temperatures  for  the  diesel  surrogate.  The  study  concluded  that 
under  HCCI  conditions,  in  which  near  stoichiometric  local  condi¬ 
tions  are  not  encountered,  any  NOx  increase  due  to  biodiesel  would 
be  small. 

3.5.  Fuel  blends  and  additives 

Apart  from  neat  fuels,  single-zone  models  have  been  used  with 
fuels  blends,  to  assess  the  effects  of  additives/blends  on  the  HCCI 
combustion  process. 

Szybist  et  al.  [38]  used  a  single-zone  model  to  simulate  the  com¬ 
bustion  of  biodiesel  blends  in  HCCI  engines.  In  their  study  they  used 
CHEMKIN  [94]  for  the  chemical  kinetics  simulation  coupled  with 
chemical  mechanisms  for  n-heptane  [99],  which  represented  diesel 
fuel,  and  methyl  butanoate  [80],  which  accounted  for  biodiesel.  The 
combined  mechanism  consisted  of  670  species  and  over  3000  reac¬ 
tions,  including  NOx  formation  reactions,  which  were  taken  from 
GRI-Mech  [84].  With  this  mechanism  it  was  possible  to  study  dif¬ 
ferent  blends  of  n-heptane  and  methyl  butanoate,  which  were  used 
as  surrogates  for  the  simulation  of  diesel/biodiesel  blends.  Thus,  the 
researchers  simulated  diesel/biodiesel  blends  of  0, 1 0  and  50%  (vol.) 
in  biodiesel  with  n-heptane/methyl  butanoate  blends,  the  relative 
proportion  of  which  was  chosen  to  match  the  oxygen  content  (wt%) 
of  the  experimental  fuel  blends. 

The  simulation  model  predicted  the  effect  of  inlet  temperature 
increase  on  the  phasing  and  magnitude  of  the  low  temperature 
heat  release  (LTHR),  and  the  phasing  of  the  main  combustion  event. 
However,  discrepancies  were  observed  between  experiments  and 


simulation  as  regards  the  dependence  of  the  results  on  the  biodiesel 
concentration,  and  the  duration  and  phasing  of  the  LTHR  and  the 
main  combustion  event.  These  were  attributed  to  the  low  cetane 
number  of  methyl  butanoate,  the  approximation  of  the  IVC  tem¬ 
perature  with  the  intake  manifold  temperature,  and  the  inherent 
inability  of  the  single-zone  model  to  capture  any  fuel  or  tempera¬ 
ture  stratification.  The  latter  resulted  to  an  overestimation  of  the 
peak  combustion  temperature  and  the  combustion  rate,  which  led 
to  almost  100%  combustion  efficiency,  very  low  CO  and  high  NOx 
emissions. 

Saisirirat  et  al.  [100]  examined  the  effect  of  n-heptane  blends 
with  ethanol  and  1 -butanol,  both  experimentally  and  numerically. 
They  used  the  SENKIN  single-zone  model  of  the  CHEMKIN  II  code 
[94]  to  simulate  n-heptane/ethanol  and  n-heptane/1 -butanol  fuel 
blends.  The  chemical  mechanism  used  was  that  of  Saisirirat  et  al. 
[44]  consisting  of  1046  species  and  4398  reactions,  which  includes 
ethanol,  1 -butanol  and  PRF  oxidation  sub-schemes.  The  simula¬ 
tion  compared  the  HCCI  combustion  of  neat  n-heptane  to  the  one 
obtained  with  the  two  fuel  blends,  i.e.  n-heptane/ethanol  and  n- 
heptane/1 -butanol  both  in  63:37  molar  proportions.  It  was  found 
that  the  addition  of  ethanol  or  1 -butanol  to  n-heptane  affected  the 
OH  radical  production  and  slowed  down  the  chemical  reactions 
during  LTHR,  decreasing  the  concentration  of  reactive  intermediate 
species  (HCHO,  OH,  H02,  etc.).  Thus,  alcohols  induced  a  chemical 
effect  on  the  LTHR.  Owing  to  this  reduction  in  chemically  active 
species  and  to  the  low  LTHR,  which  prevented  a  further  tempera¬ 
ture  rise  of  the  mixture,  the  main  combustion  event  was  delayed 
with  the  addition  of  alcohols. 

Mack  et  al.  [101]  constructed  a  single-zone  model  for  the  sim¬ 
ulation  of  ethanol/diethyl  ether  (DEE)  blends.  The  single-zone 
model  incorporated  Woschni’s  correlation  for  the  estimation  of 
heat  losses;  it  used  a  DME-derived  chemical  mechanism  [52]  for 
the  chemistry  of  DEE  oxidation,  and  Marinov’s  ethanol  oxidation 
scheme  [60]  for  the  combustion  chemistry  of  ethanol.  The  com¬ 
bined  mechanism  consisted  of  112  species  and  484  reactions.  The 
purpose  of  the  study  was  to  investigate  the  relative  reactivity  of  DEE 
and  ethanol  in  a  fuel  blend  of  25%  DEE  in  ethanol.  The  experiments 
conducted  with  carbon-14  isotope  tracing  showed  that  DEE  was 
more  reactive  than  ethanol,  yielding  more  C02  relative  to  ethanol. 
The  simulation  results  supported  this  conclusion,  since  it  was  found 
that  DEE  began  to  react  earlier  than  ethanol.  However,  a  signifi¬ 
cant  elongation  of  the  heat  release,  which  would  be  desirable  to 
control  the  combustion  process,  was  not  observed  using  the  fuel 
blend. 

The  same  research  group  also  studied  the  effect  of  small 
amounts  of  di-tertiary  butyl  peroxide  (DTBP)  (up  to  3%)  on  the  HCCI 
combustion  of  ethanol  and  diethyl  ether  blends  [45].  The  intention 
was  to  elongate  the  combustion  event  using  the  exothermic  decom¬ 
position  of  DTBP  to  promote  early  DEE  ignition.  The  single-zone 
model  used  was  the  one  mentioned  earlier  and  the  combustion 
reaction  scheme  was  enhanced  by  the  adiabatic  decomposition 
of  DTBP  based  on  the  work  of  Iizuka  et  al.  [102].  It  was  found 
via  simulation,  and  verified  by  experiments,  that  the  addition  of 
DTBP  advanced  the  ignition  timing  in  the  DEE-ethanol  blend.  The 
simulations  showed  that  the  decomposition  of  DTBP  raised  the 
temperature  of  the  mixture  and  advanced  ignition,  although  poten¬ 
tial  chemical  effects  of  DTBP  were  not  excluded. 

Yao  et  al.  [103]  used  a  zero-dimensional  thermodynamic  model 
coupled  with  CHEMKIN  III  to  simulate  a  DME/methanol  blend  in  an 
HCCI  engine.  The  chemical  mechanisms  used  were  the  ones  of  Cur¬ 
ran  et  al.  for  DME  [52],  consisting  of  78  species  and  336  reactions, 
and  Held  and  Dryer  [74]  for  the  reactions  of  methanol,  consisting  of 
22  species  and  89  reactions.  The  authors  compared  the  simulation 
predictions  to  experimental  results  of  pressure  and  operating  range 
in  %dme  -  ^MeOH  map.  The  analysis  of  the  simulation  results  showed 
that  in  DME/methanol  dual  fuel  oxidation  the  low  temperature  heat 
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release  of  DME  was  inhibited,  resulting  only  in  a  high  temperature 
heat  release.  The  OH  radicals  were  mainly  formed  by  the  decompo¬ 
sition  of  H202  at  blue-flame  reactions,  and  by  H  +  02  =  O  +  OH  and 
O  +  H20  =  OH  +  OH  at  hot  flame  reactions. 

Similar  results  were  found  by  the  same  Institute  in  another  study 
[104].  The  simulation  model  used  has  already  been  described  in  the 
first  part  of  their  paper  [48]  presented  in  Section  3.3.  In  Ref.  [104] 
the  researchers  studied  the  HCCI  combustion  of  a  DME/methane 
blend.  Under  the  simulation  conditions,  the  dual  fuel  blend  showed 
a  two  stage  autoignition  scheme.  This  was  different  from  the  clas¬ 
sical  two  stage  ignition  observed  in  neat  DME  oxidation,  since  it 
corresponded  to  two  high  temperature  heat  release  rates,  one  of 
DME,  which  was  found  to  oxidate  first,  and  the  other  of  methane. 
It  was  observed  that  methane  -  like  methanol  -  inhibited  the  LTHR 
but  promoted  the  high  temperature  heat  release  (HTHR)  of  DME. 
On  the  other  hand,  DME  promoted  methane  oxidation  both  due  to 
thermal  and  chemical  considerations.  The  thermal  effect  was  due 
to  the  temperature  increase,  which  was  the  result  of  DME  oxida¬ 
tion,  while  the  chemical  effect  was  due  to  HCHO  and  H202,  which 
were  formed  during  DME  oxidation  and  increased  the  OH  radical 
pool.  The  researchers  also  identified  the  major  reaction  paths  for 
the  dual-fuel  oxidation. 

In  [105]  Ogawa  et  al.  studied  the  suppression  effects  of  addi¬ 
tives  such  as  water,  methanol,  ethanol,  hydrogen  and  methane,  on 
the  combustion  of  a  DME  fueled  HCCI  engine.  The  researchers  used 
SENKIN  under  adiabatic  conditions  for  the  single-zone  simulations 
and  the  corresponding  chemical  mechanisms  were  adapted  from 
Curran  et  al.  [52]  for  DME,  Marinov  [60]  for  Ethanol,  and  Westbrook 
et  al.  [75,76]  for  Methanol.  In  the  experiments  the  suppressors  were 
injected  directly  into  the  combustion  chamber,  while  the  DME  was 
introduced  into  the  air  stream  at  the  intake  pipe.  The  direct  injec¬ 
tion  of  suppressors  induced  a  cooling  effect  on  the  charge.  Water 
injection  had  a  predominantly  cooling  effect,  i.e.  it  lowered  the  tem¬ 
perature  of  the  mixture  during  its  vaporization,  thereby  retarding 
ignition.  Simulations  conducted  with  the  assumption  that  gaseous 
water  was  introduced  in  the  mixture  at  BDC,  resulted  in  a  slight 
delay  in  radical  formation  and  suggested  that  the  water  addition 
also  had  a  minor  chemical  effect.  Similar  findings  characterized 
methanol  addition.  Its  main  effects  were  a  thermal  one  observed 
during  the  experiments  and  a  chemical  one  suggested  by  simula¬ 
tions. 

As  in  the  case  of  water  addition  simulations,  methanol  addi¬ 
tion  simulations  were  conducted  assuming  that  methanol  existed 
in  gaseous  state  at  the  beginning  of  calculations.  It  was  found  that 
methanol  addition  severely  suppressed  LTHR  to  near-zero  values, 
significantly  reduced  OH,  H  and  O  radicals,  and  retarded  main  com¬ 
bustion.  Further  simulations  with  ethanol,  hydrogen  and  methane 
addition  showed  that  LTHR  was  slightly  reduced  with  hydrogen  or 
methane  addition,  and  essentially  disappeared  with  the  addition 
of  ethanol  or  methanol,  in  order  of  significance.  Moreover,  it  was 
found  that  despite  the  different  simulation  conditions  with  respect 
to  additive  and  its  quantity,  the  OH  concentration  at  the  onset  of 
LTHR  remained  fairly  constant.  This  suggested  that  the  chemical 
effect  induced  by  the  additives  was  linked  to  a  reduction  in  the  OH 
radicals  prior  to  LTHR. 

This  assumption  regarding  OH  radicals  was  supported  by  later 
studies  with  H2  and  MeOH  addition  to  DME-air  mixture.  Specifi¬ 
cally,  Shudo  et  al.  [106,107]  investigated  the  effect  of  H2  addition, 
which  originated  from  on-board  methanol  reformed  gas,  on  the 
DME  HCCI  combustion.  The  concept  of  the  combustion  engine  sys¬ 
tem  and  the  experimental  setup  are  given  in  Fig.  6(a)  and  (b). 

The  SENKIN  [94]  simulation  model  was  used  to  describe  the 
effect  of  hydrogen  addition,  and  the  reaction  mechanism  used  was 
that  of  Fischer  et  al.  [53-55]  consisting  of  79  species  and  351  reac¬ 
tions.  The  simulation  results  showed  that  the  addition  of  H2  in  the 
DME-air  mixture  suppressed  DME  consumption.  This  was  due  to  an 
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Fig.  6.  (a)  Concept  of  HCCI  combustion  engine  system  fueled  with  DME  and  MRG 
onboard-reformed  from  methanol  utilizing  engine  exhaust  gas  heat  [107]  and  (b) 
experimental  setup  [106]. 


antagonistic  effect  for  the  consumption  of  OH  radicals;  in  the  pres¬ 
ence  of  H2,  OH  radicals  reacted  preferentially  with  H2  instead  of 
DME,  to  yield  the  relatively  stable  -  at  low  temperatures  -  interme¬ 
diates  H02  and  H202.  This  led  to  a  suppression  of  LTHR,  a  decrease 
in  the  in-cylinder  temperature  rise  and  a  consequent  retardation  of 
the  main  combustion  event.  This  thermal  effect  was  the  main  rea¬ 
son  for  the  combustion  retardation,  since  the  ignition  temperature 
was  not  altered  appreciably.  These  findings  implied  that  hydrogen 
or  hydrogen-containing  additives,  such  as  methanol  reformed  gas 
(MRG),  could  be  used  to  control  the  ignition  timing  of  DME.  The 
researchers  also  studied  the  conversion  of  methanol  to  DME  and  H2, 
via  dehydration-  and  thermal  decomposition-methanol  reformers, 
respectively,  to  conclude  that  the  H2/DME  fraction  in  MRG  could  be 
continuously  controlled  by  adjusting  the  methanol  supply  in  each 
reformer.  The  addition  of  CO  in  neat  DME,  which  was  studied  in 
[106],  was  also  found  to  suppress  the  DME  oxidation  at  low  tem¬ 
perature,  but  its  effect  was  less  pronounced  than  the  H2  addition 
effect. 

Yamada  et  al.  [108]  studied  the  effect  of  methanol  or  03  addi¬ 
tion  to  DME-air  mixture  under  HCCI  combustion  conditions.  They 
used  the  SENKIN  [94]  code  to  simulate  the  engine  operation,  and 
implemented  a  chemical  mechanism  by  Curran  et  al.  [52  ]  to  account 
for  the  DME  and  methanol  addition  chemistry.  Ozone  addition  was 
accounted  for  by  including  the  thermal  decomposition  reaction  of 
ozone  according  to  Heirmel  and  Coffee  [109].  The  experimental 
and  simulation  results  showed  that  methanol  addition  retarded  the 
main  heat  release  event  significantly;  8%  addition  of  methanol  in 
DME  resulted  in  a  1 5°  CA  delay.  It  was  suggested  that  this  delay  was 
due  to  the  competitive  consumption  of  OH  by  methanol  instead  of 
DME,  which  led  to  the  formation  of  relatively  stable  intermediates 
such  as  HCHO  and  H02.  In  contrast,  the  addition  of  low  03  amounts 
(0.01 5%  in  DME)  advanced  the  ignition  timing  by  as  much  as  20°  CA. 
This  was  due  to  the  decomposition  of  03  at  relatively  low  temper¬ 
atures  (~600I<),  which  provided  oxygen  atoms  and  activated  the 
mixture.  This  supply  of  oxygen  atoms  at  temperatures  lower  than 
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those  of  normal  cool  flames  suppressed  the  production  of  formalde¬ 
hyde,  which  was  a  chain  termination  step  under  low  temperatures. 
Thus,  more  heat  was  liberated  during  cool  flame  reactions,  leading 
to  an  earlier  thermal  ignition. 

4.  Multi- zone  modeling 

The  limitations  encountered  with  single-zone  models  can  be 
partly  overcome  with  the  use  of  multi-zone  models.  The  latter 
assume  that  the  combustion  chamber  is  divided  into  a  sufficient 
number  of  regions  or  zones,  each  possessing  uniform  thermody¬ 
namic  properties.  The  sufficiency  of  the  number  of  zones  is  usually 
determined  by  a  sensitivity  analysis,  i.e.  the  minimum  number  of 
zones  used  is  defined  as  the  number  beyond  which  the  simulation 
results  tend  to  converge. 

The  temperature,  composition,  volume  and  all  thermodynamic 
properties  of  each  zone  are  considered  uniform  within  the  bound¬ 
aries  of  the  zone,  but  can  be  different  from  the  thermodynamic 
properties  of  other  zones.  This  provides  an  opportunity  to  study 
the  effects  of  temperature  or  species  stratification  within  the  com¬ 
bustion  chamber,  and  usually  leads  to  a  better  estimation  of  the 
combustion  duration,  the  pressure  rise  rate  and  the  emissions  for¬ 
mation.  Mass  and  heat  transfer  within  the  combustion  chamber 
play  a  critical  role  in  the  estimation  of  the  aforementioned  operat¬ 
ing  characteristics  and  since  the  combustion  chamber  is  considered 
an  ensemble  of  thermodynamically  distinct  zones,  the  processes 
of  mass  and  heat  transfer  between  zones  must  also  be  described 
by  the  multi-zone  model.  Fig.  7  shows  the  first  law  of  thermo¬ 
dynamics  applied  to  a  zone,  which  in  the  general  case  includes 
heat  and  mass  transfer.  As  regards  mass  transfer,  each  zone  can  be 
thermodynamically  defined  as  a  closed  or  an  open  thermodynamic 
system,  depending  on  whether  mass  flow  between  zones  is  pro¬ 
hibited  or  allowed,  respectively.  Likewise,  heat  transfer  between 
zones  and  to  the  combustion  chamber  wall  can  also  be  permitted 
or  neglected. 

However,  the  assumptions  regarding  mass  and  heat  transfer 
affect  the  results  of  the  simulation  significantly;  heat  transfer 
affects  the  temperature  stratification  within  the  combustion  cham¬ 
ber,  the  colder  regions  of  which  are  sources  of  unburned  HC  and  CO. 
The  formation  of  these  pollutants  is  also  affected  by  mass  trans¬ 
fer  between  zones,  since  any  HC  or  CO  flowing  from  a  colder  to  a 
hotter  region/zone  is  expected  to  be  (partially)  oxidized  [110].  The 
predicted  combustion  rate,  pressure  rise  rate,  combustion  duration 
and  peak  pressure  are  also  affected  by  the  choice  of  the  sub¬ 
models  used  to  describe  heat  and  mass  transfer.  The  multi-zone 
models  examined  subsequently  are  divided  into  stochastic  based 
multi-zone  models,  in  which  some  properties  of  the  mixture  are 
considered  random  variables,  and  phenomenological  multi-zone 
models,  in  which  probability  considerations  have  been  excluded. 


In  [62]  Mosbach  et  al.  constructed  a  stochastic  reactor  model 
(SRM)  based  on  a  probability  density  function  (PDF)  approach.  The 
PDF  was  approximated  by  an  ensemble  of  100  stochastic  particles. 
These  did  not  contain  any  spatial  or  geometric  information.  The 
random  variables  were  the  species  concentration  and  the  tempera¬ 
ture.  The  model  included  heat  transfer  to  the  combustion  chamber 
wall  using  Woschni’s  correlation  [91]  and  the  magnitude  of  the 
temperature  fluctuations  was  regulated  by  a  coefficient.  Turbulent 
mixing  was  accounted  for  via  a  mixing  intensity  term  that  included 
the  characteristic  mixing  time.  The  fuel  simulated  was  a  blend  of 
EtOH/DEE,  although  a  validation  of  the  model  was  also  presented 
with  neat  EtOH  as  fuel.  The  validation  showed  good  agreement 
between  experimental  and  simulated  results  as  regards  pressure 
traces.  The  model  also  captured  qualitatively  the  experimental  vari¬ 
ation  of  CO  and  HC  emissions  with  the  combustion  phasing,  i.e. 
the  angle  at  which  50%  of  the  total  heat  has  been  released.  In 
these  simulations,  the  equivalence  ratio  remained  constant  and 
the  combustion  phasing  was  regulated  via  the  initial  temperature. 
The  stochastic  reactor  model  was  also  used  to  study  numerically 
the  effect  of  the  DEE  in  EtOH  ratio  on  the  combustion  phasing, 
the  combustion  duration  and  the  emissions  of  the  HCCI  engine.  It 
was  shown  that  as  the  fraction  of  DEE  in  the  DEE/EtOH  blend  was 
increased,  the  combustion  phasing  advanced  before  TDC,  the  com¬ 
bustion  duration  was  shorter,  and  the  CO  emissions  were  reduced. 
This  fact  provided  indication  that  DEE  could  be  used  with  ethanol 
to  regulate  the  combustion  phasing  and  duration  of  the  HCCI  com¬ 
bustion. 

Noda  and  Foster  [1 1 1  ]  numerically  studied  hydrogen  HCCI  com¬ 
bustion  by  constructing  a  multi-zone  model;  the  different  zones 
catered  for  any  temperature  or  composition  inhomogeneities.  The 
zones  were  not  assigned  a  specific  location  within  the  combus¬ 
tion  chamber,  and  mass  exchange  between  them  was  neglected. 
The  heat  losses  were  calculated  via  the  Woschni  [91]  correlation 
using  the  combustion  chamber  area  and  the  mean  gas  tempera¬ 
ture.  The  combustion  chemistry  was  accounted  for  by  incorporating 
11  species  and  23  reactions  from  Westbrook  and  Dryer  [69]  for 
the  oxidation  of  hydrogen,  and  the  extended  Zeldovich  mecha¬ 
nism  of  Lavoie  et  al.  [83]  for  NOx  formation.  The  researchers  used 
10-30  zones  to  examine  the  effect  of  temperature  and  fuel  inho¬ 
mogeneities  on  HCCI  combustion.  Both  temperature  and  -  to  a 
much  lesser  degree  -  fuel  inhomogeneities  increased  the  com¬ 
bustion  duration  compared  to  single-zone  simulations.  The  mean 
gas  temperature  and  the  type  of  temperature  inhomogeneity  also 
affected  combustion  duration  and  the  formation  of  NO  emissions. 
Provided  the  mean  gas  temperature  was  high  enough  to  avoid 
bulk  quenching  of  the  mixture,  the  most  favorable  type  of  tem¬ 
perature  stratification  included  a  small  amount  of  hot  gases  to 
induce  autoignition  and  avoid  excessive  NO  formation  and  a  greater 
amount  of  colder  mixture.  This  type  of  stratification  expanded  the 
operating  range  in  a  mean  gas  temperature  versus  temperature 
inhomogeneity  map. 

Komninos  et  al.  [112]  and  Komninos  and  Rakopoulos  [40,113] 
used  a  multi-zone  model  to  simulate  hydrogen  [112]  and  ethanol 
[40,113]  HCCI  combustion.  In  all  studies,  the  multi-zone  model 
used  divided  the  combustion  chamber  into  several  zones,  which 
occupied  specific  geometric  regions  (Fig.  8). 

The  only  uniformity  assumed  was  that  of  pressure,  while  heat 
and  mass  transfer  between  zones  was  accounted  for,  taking  into 
account  the  geometric  configuration  of  the  zones;  each  zone  could 
communicate  with  its  immediate  inner  and  outer  zone  (Fig.  9).  The 
model  also  accounted  for  a  crevice  zone,  which  was  in  direct  com¬ 
munication  with  the  outmost  zone  of  the  combustion  chamber.  In 
Ref.  [  1 1 2  ],  1 1  zones  were  used  -  including  the  crevices  -  to  simulate 
hydrogen  HCCI  combustion  and  the  Annand  correlation  [114]  was 
employed  to  estimate  the  heat  flux  from  the  outmost  zone  to  the 
combustion  chamber  wall.  The  chemistry  of  hydrogen  oxidation 
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Fig.  8.  Zone  configuration  in  the  multi-zone  model  [112]. 


was  accounted  for  via  the  mechanism  of  Marinov  et  al.  [54]  involv¬ 
ing  10  species  and  19  reactions,  augmented  by  a  NOx  formation 
scheme  adapted  from  the  isooctane  oxidation  reaction  mechanism 
of  Golovitchev  [86].  A  parametric  investigation  was  conducted  with 
the  multi-zone  model  to  determine  the  effect  of  various  operating 
parameters  on  hydrogen  HCCI  combustion.  Compression  ratio  and 
initial  temperature  were  the  dominant  factors  affecting  the  igni¬ 
tion  timing  and  pressure  rise  rate.  Engine  speed  also  affected  the 
ignition  timing,  mainly  due  to  the  variation  of  the  available  time  for 
the  occurrence  of  chemical  reactions.  It  was  found  that  increasing 
the  load,  i.e.  the  equivalence  ratio,  NOx  emissions  were  increased 
significantly.  Inlet  pressure  increased  almost  proportional  to  the 
attainable  imep,  without  compromising  the  fuel  conversion  effi¬ 
ciency  for  the  conditions  simulated. 

In  [40]  Komninos  and  Rakopoulos  numerically  compared 
ethanol  and  isooctane  HCCI  combustion.  The  multi-zone  model 
used  was  the  one  described  in  the  preceding  paragraph  with  some 
modifications;  the  zones  number  was  increased  to  16,  the  zone 
resolution  was  refined  near  the  combustion  chamber  wall,  and  the 
wall  heat  flux  estimation  approach  abandoned  the  Annand  corre¬ 
lation  and  was  based  on  the  temperature  gradient  near  the  wall. 
These  modifications  were  necessary  to  better  estimate  the  forma¬ 
tion  of  HC  and  CO  emissions  [115].  The  chemical  mechanism  used 
for  the  ethanol  oxidation  was  that  of  Rohl  and  Peters  [57]  consist¬ 
ing  of  38  species  and  228  reactions.  Since  this  mechanism  did  not 
include  NOx  reactions,  it  was  augmented  by  13  reactions  involving 
nitrogen  oxides  formations  [86].  The  model  vas  validated  against 


experimental  results  and  captured  adequately  the  pressure  traces 
and  HRR  for  both  ethanol  and  isooctane. 

According  to  the  simulation  results,  CO  emissions  were  formed 
near  the  combustion  chamber  walls,  thereby  justifying  the  finer 
zonal  resolution  used  in  these  regions  (see  Fig.  10).  As  shown  in 
Fig.  10,  the  amount  and  location  of  CO  production  was  found  to 
depend  on  load;  at  higher  loads,  CO  production  advanced  to  the 
inner  zones. 

Total  CO  emissions  were  adequately  captured  but  HC  emissions 
were  overestimated,  especially  for  ethanol.  An  investigation  for  the 
composition  of  the  HC  emissions  for  both  isooctane-  and  ethanol- 
fueled  HCCI  combustion  showed  that  in  the  ethanol  cases  examined 
the  vast  majority  of  unburned  HC  consisted  of  oxygenated  HC 
species  (OHC).  These  were  primarily  ethanol,  acetaldehyde  and 
formaldehyde.  The  predicted  percentage  of  OHC  in  total  HC  (aprox. 
95%  by  mass)  for  ethanol  was  not  very  sensitive  to  the  load 
variation  under  the  simulation  conditions.  On  the  contrary,  the 
relative  amount  of  OHC  in  total  HC  was  lower  (19-33%  by  mass) 
in  the  case  of  isooctane  and  decreased  with  the  increase  of  load. 
It  was  concluded  and  demonstrated  that  this  would  induce  a 
significant  error  in  the  HC  measurement  if  a  measuring  device 
such  as  the  FID  was  used,  since  the  FID  principle  possesses  an 
inherent  inability  to  accurately  detect  oxygenated  carbon  atoms 
[116]. 

Using  the  same  model  in  [113]  the  authors  compared  two 
ethanol  oxidation  mechanisms,  namely  Marinov’s  [60]  and  Rohl 
and  Peters’  [57].  Marinov’s  mechanism  consisted  of  57  species  and 
383  reactions,  while  Rohl  and  Peters’  mechanism,  which  was  a 
reduced  version  of  Marinov’s,  consisted  of  38  species  and  228  reac¬ 
tions.  Both  mechanisms  were  augmented  by  NOx  related  reactions 
taken  from  [86].  The  results  showed  that  the  two  mechanisms 
predicted  essentially  the  same  dominant  HC  species,  with  the 
exception  of  formic  acid  (HCOOH)  and  to  a  lesser  extent  methanol 
(CH3OH),  both  of  which  were  not  included  in  the  reduced  mech¬ 
anism.  The  mechanism  reduction  had  a  significant  impact  on 
the  computational  time:  although  the  species  used  in  Rohl  and 
Peters’  mechanism  were  2/3  of  the  species  used  in  the  Marinov 
mechanism,  the  Rohl  and  Peters’  mechanism  required  half  the  com¬ 
putational  time.  The  simulation  results  showed  a  slightly  higher 
predicted  combustion  rate  by  the  Rohl  and  Peters’  mechanism  rel¬ 
ative  to  Marinov’s  mechanism,  despite  the  fact  that  the  predicted 
initial  heat  release  was  the  same  for  the  two  mechanisms.  The 
resulting  temperature  fields  were  almost  identical  as  shown  in 
Fig.  11.  Therefore,  the  two  chemical  mechanisms  did  not  induce 
significant  differences  in  the  predicted  performance.  NOx,  CO  and 
HC  emissions  were  essentially  the  same  for  the  two  mechanisms, 
with  the  greatest  disagreement  occurring  between  predicted  CO 
emissions  (~9%)  at  the  low  load  cases.  This  was  connected  with  the 
slightly  lower  combustion  rate  predicted  by  Marinov. 

Sjoberg  and  Dec  [117,118]  used  a  multi-zone  model  to  study  the 
effects  of  engine  speed,  equivalence  ratio,  initial  temperature  and 
boost  pressure  on  ethanol  HCCI  combustion.  They  used  the  multi¬ 
zone  model  provided  by  the  Senkin  application  of  the  CHEMKIN 
package  [94],  which  treated  each  zone  as  a  single  lumped  mass 
with  uniform  composition.  The  only  interaction  between  zones 
was  through  pressure,  which  was  considered  uniform.  Thus,  the 
volume  of  each  zone  depended  on  its  mass  and  its  temperature 
relative  to  the  other  zones.  The  1 1  zones  used  were  adiabatic  in 
principle.  This  led  to  the  need  to  introduce  initial  temperature 
stratification  for  the  10  active  zones  to  match  the  experimental 
pressure  trace.  The  11th  zone  was  assigned  a  temperature  low 
enough  to  render  it  chemically  inactive  throughout  the  engine 
cycle.  This  zone  included  16%  of  the  total  mass  and  compensated 
for  the  combustion  inefficiencies,  heat  transfer  and  blow-by,  which 
were  not  directly  modeled.  The  simulated  compression  ratio  was 
also  reduced  relative  to  the  experimental  one  to  account  for  the 
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Fig.  10.  CO  simulation  results  for  ethanol;  left:  in-cylinder  CO  emissions  and  right:  zone  CO  production  rate  [40]. 


neglected  heat  losses.  The  thermal  stratification  of  the  zones  was 
supported  by  KIVA  simulations  and  reproduced  adequately  the 
variation  of  combustion  phasing  with  temperature.  The  chemical 
mechanism  used  consisted  of  58  species  and  310  reactions  [61], 
and  was  part  of  a  larger  gasoline  surrogate  mechanism  produced 
by  the  joint  effort  of  the  National  University  of  Ireland-Galway 
and  Lawrence  Livermore  National  Laboratories  [119]. 

The  simulations  above  showed  that  the  ethanol  oxidation  mech¬ 
anism  reproduced  adequately  the  effect  of  equivalence  ratio,  engine 
speed  and  boost  pressure  on  ignition.  The  same  model  was  used  by 
these  researchers  to  study  the  effect  of  EGR  and  its  constituents  on 
the  HCCI  autoignition  of  ethanol  [118].  In  this  study,  it  was  found 
that  the  addition  of  simulated  EGR  did  not  affect  significantly  the 
ethanol  autoignition  timing  for  the  conditions  examined.  The  effect 
of  EGR  on  other  fuels,  such  as  isooctane  and  PRF,  was  more  pro¬ 
nounced.  This  was  mainly  owed  to  the  low  sensitivity  of  ethanol 
oxidation  on  oxygen  concentration,  which  was  reduced  with  the 


increase  of  simulated  EGR.  Moreover,  it  was  suggested  that  the 
presence  of  water  in  the  EGR  promoted  ethanol  auto-ignition. 

5.  Multi-dimensional  modeling 

Besides  the  aforementioned  favorable  attributes  of  the  multi¬ 
zone  models,  they  also  have  limitations.  Even  in  the  cases  where 
spatial  considerations  have  been  accounted  for  by  multi-zone  mod¬ 
els,  they  usually  do  not  allow  high  local  resolution.  Moreover, 
there  is  an  inherent  difficulty  to  fundamentally  describe  the  motion 
of  the  charge  within  the  combustion  chamber  and  the  effects  of 
turbulence.  The  problem  is  even  more  pronounced  in  the  case 
of  Diesel-like  HCCI  engines,  in  which  the  mixture  formation  is 
achieved  via  direct  fuel  injection  within  the  combustion  cham¬ 
ber.  Direct  fuel  injection  imposes  new  demands  on  any  simulation 
model,  since  the  mixture  -  and  the  resulting  temperature  -  strati¬ 
fication  within  the  combustion  chamber  is  affected  significantly  by 
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Fig.  11.  Zone  temperature  distributions  produced  by  Rohl  and  Peters’  (solid  line) 
and  Marinov’s  (dashed  line)  mechanisms  for  the  high  load  ethanol  case  (0  =  0.2641 ) 
[113]. 

the  injection  characteristics.  Phenomena  such  as  wall  wetting,  fuel 
dispersion,  fuel  vaporization  under  the  low  pressure-low  temper¬ 
ature  conditions  encountered  during  early  direct  injection,  must 
be  described.  In  these  cases,  the  incorporation  of  accurate  models 
for  the  estimation  of  the  physicochemical  fuel  properties  becomes 
crucial,  especially  for  biofuels,  since  their  development  and  use  has 
a  relatively  short  history. 

The  aforementioned  issues  can  be  dealt  with  by  using  multi¬ 
dimensional  models,  i.e.  CFD  models,  which  describe  in  detail  the 
charge  motion  within  the  combustion  chamber  and  the  turbulence 
effects,  provide  a  superior  spatial  resolution  of  the  combustion 
chamber,  and  can  account  for  the  description  of  DI HCCI  combustion 
by  incorporating  the  dynamics  and  kinematics  of  the  fuel  injection. 

In  the  remaining  part  of  this  section,  direct  injection  HCCI  CFD 
models  will  be  discussed  first  while  dedicated  HCCI  CFD  models 
assuming  homogeneous  conditions  at  IVC  will  be  presented  subse¬ 
quently. 

5.2.  Direct  injection  HCCI 

Kim  et  al.  [120]  investigated  the  influence  of  the  spray  angle 
and  advanced  injection  timing  on  the  combustion  performance 
and  emission  characteristics  in  a  DME  fueled  compression  ignition 
engine.  They  used  the  KIVA  software  [121],  which  was  supple¬ 
mented  with  liquefied  DME  fuel  properties  and  their  temperature 
dependence  [122-126].  The  oxidation  process  of  the  DME  fuel 
was  described  with  a  detailed  chemical  kinetic  model  [53-55], 
which  consisted  of  351  reversible  elementary  reactions  among  79 
species.  This  mechanism  was  expanded  by  including  NOx  forma¬ 
tion  reactions  taken  from  a  modified  version  of  GRI  Mech  [84].  A 
two-step  phenomenological  soot  model  [127]  was  also  applied  to 
calculate  soot  emissions.  The  DME  spray  and  atomization  char¬ 
acteristics  were  simulated  via  the  Kelvin-Helmholtz  (KH)  and 
Rayleigh-Taylor  (RT)  hybrid  breakup  model  [128].  Turbulence  was 
accounted  for  via  the  modified  Renormalization  Group  (RNG)  k-s 
model  [129],  and  wall  impingement  was  simulated  via  the  spray 
wall  interaction  model  [130]  including  droplet  splash/spreading.  A 
computational  grid  with  a  unit  cell  size  of  about  2  mm3  was  formed 
for  the  evaluation  of  spatial  distribution  in  the  cylinder  and  the 
computational  grid  was  divided  into  a  60°  sector  mesh,  due  to  the 
symmetric  arrangement  of  the  six-hole  injector.  The  calculation 
time  step  was  set  at  0.1  |jis. 


Fig.  12.  Cylinder  locations  according  to  the  injection  timing  and  five  injection  angles 
[120]. 

The  model  was  validated  against  experimental  measurements 
at  three  injection  timings,  with  an  injection  angle  of  156°.  The 
experimental  and  simulation  results  at  an  injection  timing  of  60° 
BTDC  had  two  peaks  in  the  heat  release  rate,  indicating  low-  and 
high-temperature  DME  oxidation.  The  validated  model  was  subse¬ 
quently  used  for  a  numerical  study  aimed  at  analyzing  the  effect 
of  the  injection  angles  and  timings  shown  in  Fig.  12.  The  result¬ 
ing  fuel  spray  patterns  created  by  these  configurations  are  shown 
in  Fig.  13.  These  patterns  revealed  fuel  spray- wall  interactions  and 
accumulation  of  fuel  in  the  squish  area  and,  possibly,  the  crevice 
regions. 

Generally,  the  fuel  injected  vaporized  rapidly  as  it  entered  the 
combustion  chamber  due  to  the  high  volatility  of  DME.  At  an  early 
injection  timing  (60°  bTDC),  the  variation  of  the  injection  angle 
affected  significantly  the  pressure  and  heat  release  rate  diagrams; 
at  wider  injection  angles  the  peak  pressure  and  HRR  were  both 
reduced,  since  some  of  the  injected  fuel  reached  the  relatively  cold 
crevice  and  combustion  chamber  wall  regions.  At  later  injection 
timings,  the  effect  of  the  injection  angle  was  much  less  pronounced, 
since  the  fuel  was  injected  within  the  piston  bowl.  As  the  injection 
timing  approached  TDC,  the  combustion  mode  shifted  progres¬ 
sively  from  partially  homogeneous  to  diffusive  combustion.  Soot 
emissions  were  very  low  for  all  cases  examined  due  to  the  chemi¬ 
cal  structure  of  the  fuel,  which  excludes  C-C  bonds.  Another  reason 
was  the  relatively  premixed  nature  of  the  mixture,  especially  at 
early  injection  timings. 

NO*  emissions  were  very  low  for  the  early  injection  timings  (IT 
40°  bTDC  and  earlier)  due  to  the  premixed  nature  of  the  mixture. 
As  the  injection  timing  approached  TDC,  NOx  emissions  increased, 
probably  owing  to  the  predominance  of  the  in-cylinder  inhomo¬ 
geneities.  At  advanced  injection  timings  and/or  wide  injection 
angles,  the  HC  and  CO  emissions  were  higher  due  to  the  introduc¬ 
tion  of  fuel  into  the  squish  and  crevice  regions.  At  narrow  angles,  the 
high  HC  and  CO  emissions  were  due  to  the  fuel  impingement  on  the 
piston  bowl.  As  regards  performance,  a  narrow  angle  produced  low 
ISFC  at  relatively  advanced  injection  timing.  Other  studies  focusing 
solely  on  the  spray  characteristics  of  DME  by  the  research  group  of 
the  Hanyang  University  can  be  found  in  [131,132] 

Chen  et  al.  [133]  investigated  experimentally  and  numerically 
the  effect  of  methanol  injection  timing  in  a  DME/MeOH  dual 
fuel  engine,  in  which  the  DME  was  introduced  at  the  inlet  port. 
For  the  numerical  simulation  of  the  closed  part  of  the  engine 
operation,  the  researchers  coupled  chemical  kinetics  with  the  CFD 
code  FIRE  [134].  The  chemical  mechanism  used  to  estimate  the 
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Fig.  13.  Effect  of  the  injection  angle  and  timing  on  calculated  spray  and  evaporation 
distributions  in  the  cylinder  according  to  the  crank  angle  [120]. 


methanol/DME  reaction  rates  consisted  of  30  species  and  44  reac¬ 
tions  including  the  fuel  chemistry  [50],  the  Zeldovich  mechanism, 
and  the  N20  intermediate  mechanism  [82].  The  effect  of  turbulence 
was  accounted  for  by  adopting  a  compressible  version  of  the  stan¬ 
dard  two-equation  k-s  turbulence  model.  The  source  terms  in  the 
transport  equations  of  FIRE  were  calculated  with  the  assumption 
that  the  reaction  rate  is  mainly  determined  by  kinetic  and  turbulent 
time  scales.  Major  submodels  used  in  the  computation  included 
the  WALLJET1  wall  interaction  model,  the  Dukowicz  evaporation 
model,  and  the  HuhGosman  breakup  model  [135].  The  Han-Reitz 
wall  heat  transfer  model  was  used  to  estimate  the  wall  heat  flux 
[136].  The  spray  angle  was  152°  formed  by  a  7-hole  injector,  which 
allowed  for  the  simulation  of  a  51°  sector  consisting  of  30,000 
grid  cells  at  BDC.  The  simulations  examined  the  effect  of  MeOH 
injection  at  two  injection  timings,  namely  -26°  and  -6°  aTDC. 

It  was  found  that  methanol  injected  at  -26°  aTDC  had  no 
effect  on  the  low  temperature  reaction  of  DME  but  accelerated  the 
heat  release  of  the  high  temperature  reaction,  producing  locally 
high  temperature  regions  and  increasing  NO  concentration.  In  this 
case,  the  injected  methanol  was  completely  oxidized.  At  the  later 
injection  timing  (-6°  aTDC)  the  DME  found  in  the  crevices  was 
higher,  and  the  high  temperature  region  was  dispersed  in  the 
clearance  volume  and  the  region  near  the  combustion  chamber 
wall,  leading  to  lower  NO  concentration.  Some  methanol  vapor 
remained  in  the  combustion  chamber  prolonging  the  combustion 
duration.  This  led  to  an  investigation  aimed  at  decreasing  the  com¬ 
bustion  duration  at  this  -  relatively  late  -  injection  timing.  The 
researchers  studied  numerically  the  effect  of  injection  pressure  on 
the  combustion  duration,  while  maintaining  the  injection  timing 


Fig.  14.  Computational  mesh  of  combustion  chamber  at  TDC  used  in  [138]. 


at  -6°  aTDC.  They  found  that  increasing  the  injection  pressure  from 
20  to  30  MPa  shortened  the  combustion  duration,  increased  the 
thermal  efficiency,  and  induced  an  increase  in  peak  temperature 
and  NO  concentration.  However,  it  was  difficult  to  experimentally 
examine  the  effect  of  increased  injection  pressure  due  to  the  low 
viscosity  of  methanol. 

5.2.  Fully  premixed  HCCI 

In  Ref.  [137]  Yao  et  al.  examined  the  HCCI  combustion  of  a  dual 
fuel  engine  fueled  with  DME/MeOH  blends.  They  used  the  STAR- 
CD/KINETICS  CFD  code  to  simulate  the  combustion  of  the  charge 
that  was  uniform  at  IVC.  The  RNG  k-s  model  was  included  to 
account  for  the  effects  of  turbulence,  and  the  PISO  algorithm  was 
used  for  the  transient  flow  of  the  engine.  The  CFD  code  was  coupled 
with  the  CHEMKIN  code  to  obtain  the  species  concentration  evolu¬ 
tion,  and  the  DME/methanol  oxidation  was  described  by  a  reduced 
mechanism  consisting  of  27  species  and  35  reactions  [50].  The  time 
step  corresponded  to  0.1°  CA  at  an  engine  speed  of  1400  rpm.  The 
researchers  validated  the  model  results  against  an  experimentally 
obtained  pressure  diagram  and  found  adequate  agreement  with 
the  measurements.  Under  the  same  experimental  conditions  they 
also  conducted  an  investigation  into  the  formation  of  HC  and  CO 
emissions.  They  found  that  UHC  emissions  mainly  consisted  of  the 
unburned  fuels  and  CH20.  The  unburned  DME  and  methanol  mainly 
resided  in  the  piston-ring  crevice  region,  while  most  of  the  CH20 
emissions  were  located  next  to  the  cylinder-liner  wall.  CO  emis¬ 
sions  were  mainly  found  in  the  regions  near  the  top  surface  of  the 
piston. 

The  authors  also  studied  the  effect  of  load  (in  terms  of  equiva¬ 
lence  ratio)  and  the  relative  proportion  of  DME  in  the  fuel  blend  on 
the  HC  and  CO  emissions.  Increasing  the  DME  ratio,  while  main¬ 
taining  the  overall  equivalence  ratio,  improved  the  combustion 
efficiency  and  reduced  the  HC  and  CO  emissions.  Moreover  at  all 
loads  examined,  a  critical  value  for  the  DME  ratio  was  found;  when 
reducing  the  ratio  below  this  value  induced  a  significant  increase  in 
both  CO  and  HC  due  to  incomplete  combustion.  However,  the  criti¬ 
cal  DME  ratio  was  not  the  same  for  all  equivalence  ratios  examined. 
At  high  equivalence  ratios,  the  critical  DME  ratio  shifted  toward 
higher  values  due  to  the  increase  in  methanol  concentration  at 
high  loads,  which  induced  an  inhibiting  effect  on  DME  ignition. 
A  critical  value  was  also  found  for  the  peak  combustion  temper¬ 
ature  (1400  K),  above  which  the  HC  and  CO  emissions  were  greatly 
reduced.  When  the  load,  i.e.  the  overall  equivalence  ratio,  was 
increased  while  maintaining  the  DME  proportion,  HC  and  CO  emis¬ 
sions  were  reduced. 

Using  roughly  the  same  simulation  model,  Huang  et  al.  [138] 
studied  numerically  the  HCCI  combustion  of  a  DME  fueled  engine. 
A  reduced  chemical  mechanism  for  DME  HCCI  combustion  pro¬ 
cess  was  adopted,  which  was  proposed  by  Yao  et  al.  [48]  and  Liang 
et  al.  [49].  The  reduced  chemical  mechanism  was  derived  from  the 
detailed  chemical  mechanism  and  consisted  of  26  species  and  28 
reactions.  A  single  layer  grid  was  applied  with  an  average  cell  size 
of  about  2  mm,  yielding  about  1658  cells  at  IVC  in  the  combustion 
chamber  (Fig.  14). 
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Fig.  15.  Initial  locations  for  LTR  and  HTR  in  [138]. 


The  computational  time  for  the  closed  part  of  the  engine  cycle 
was  approximately  18  h  in  a  Pentium  IV  3.0  GHz  CPU  PC  with  1  GB 
memory.  Under  homogeneous  charge  IVC  conditions,  the  DME  HCCI 
combustion  process  does  not  take  place  in  the  entire  combustion 
chamber  simultaneously.  The  low  temperature  reaction  was  initi¬ 
ated  primarily  at  the  region  near  the  piston  surface  and  the  squish 
region,  and  subsequently  moved  toward  the  cylinder  head.  The 
high  temperature  reaction  was  initiated  in  the  combustion  cham¬ 
ber  core  zone  and  squish  region,  and  next  proceeded  to  the  rest  of 
the  combustion  chamber  (Fig.  15). 

The  majority  of  HC  emission  was  unburned  fuel  and  CH20, 
and  resided  with  CO  in  the  bottom,  middle,  and  upper  part  of 
the  piston-ring  crevice  regions,  respectively.  At  very  low  loads, 
DME  was  partially  oxidized  producing  unburned  HC  and  CO.  As 
the  load  increased,  the  unburned  HC  decreased  and  the  CO  emis¬ 
sions  increased  and  then  decreased.  This  behavior  can  be  explained 
considering  the  successive  steps  through  which  the  hydrocarbon 
oxidation  takes  place,  i.e.  the  oxidation  follows  the  general  scheme 
HC  ->  CO  ->  C02.  At  extremely  low  loads  the  fuel  and  any  unburned 
HC  produced  by  its  decomposition  fails  to  be  converted  to  CO  to  a 
great  extent.  At  higher  loads,  the  HC  to  CO  conversion  takes  place 
more  vividly  and,  at  the  highest  attainable  loads,  both  HC  and  CO  are 
oxidized  to  produce  C02  through  the  aforementioned  generalized 
oxidation  scheme. 

Chun-Ian  et  al.  [56]  simulated  neat  DME  HCCI  combustion.  They 
used  a  single-zone  model  to  compare  two  chemical  mechanisms: 
a  detailed  mechanism  of  97  species  and  457  reactions,  which  was 
obtained  by  combining  the  mechanism  from  Curran  et  al.  [54]  with 
the  NOx  formation  reactions  of  Marinov  [85],  and  a  simplified  reac¬ 
tion  mechanism  of  36  species  and  73  reactions.  The  latter  included 
a  low  temperature,  a  negative  temperature  coefficient,  and  a  pyrol¬ 
ysis  and  oxidation  submodel  along  with  the  reaction  path  for  NOx 
formation. 

The  FLUENT  CFD  software  was  also  used  to  obtain  insight  into 
the  formaldehyde  and  formic  acid  formation  process  during  the 
DME  HCCI  combustion.  For  this  reason  the  CFD  code  was  coupled 
with  the  three  submodels,  i.e.  the  low  temperature,  the  negative 
temperature  coefficient,  and  the  pyrolysis  and  oxidation  submod¬ 
els.  The  RNG  k-e  turbulence  model  and  the  EDC  combustion 
model  were  implemented  in  the  CFD  code.  Using  this  model,  the 
researchers  found  that  both  formaldehyde  and  formic  acid  existed 
in  the  exhaust.  Their  concentration  was  close  to  the  experimen¬ 
tally  determined,  in  the  case  of  formaldehyde,  but  the  measured 
formic  acid  was  much  less  than  the  one  predicted.  The  CFD  cal¬ 
culations  showed  that  the  formation  of  these  species  was  located 
at  the  region  near  the  combustion  chamber  wall  and  above  the 
ring  crevices.  Their  formation  originated  from  the  partial  oxida¬ 
tion  of  DME,  which  had  entered  the  crevices  during  combustion 
and  flowed  back  to  the  combustion  chamber  during  expansion.  The 
partial  oxidation  of  DME  gave  rise  to  formaldehyde  and  formic  acid, 
the  further  oxidation  of  which  was  suppressed  due  to  the  low  local 
temperatures  induced  by  the  gas  expansion  and  heat  transfer  to  the 
cylinder  wall.  Some  of  the  formaldehyde  and  formic  acid  formed 
were  transferred  to  the  central  part  of  the  combustion  chamber, 


due  to  mass  transfer,  indicating  that  their  concentration  can  be 
affected  by  the  combustion  chamber  configuration  and  the  exist¬ 
ing  flow  field.  At  EVO  a  relatively  high  concentration  of  these  gases 
remained  near  the  wall. 

In  Ref.  [139]  Park  numerically  analyzed  the  combustion  and 
emission  characteristics  and  the  optimal  operating  conditions  of 
homogeneous  charge  compression  ignition  (HCCI)  engines  fueled 
with  dimethyl  ether  (DME).  For  this  investigation,  a  modified 
version  of  KIVA-3V  code  [121]  was  used  to  simulate  DME  HCCI 
combustion  with  a  two-dimensional  computational  mesh.  The 
CHEMKIN  chemistry  solver  [70]  was  integrated  into  KIVA-3V,  in 
order  to  solve  the  chemistry  and  the  chemical  mechanism  used  for 
the  simulation  of  DME  consisted  of  351  reactions  and  79  species 
[53-55].  A  two-step  phenomenological  model  was  used  to  simu¬ 
late  soot  emission,  since  the  study  extended  the  investigation  into 
rich  DME  combustion.  The  soot  model  considered  both  soot  for¬ 
mation,  using  a  modified  Hiroyasu  model  that  considers  acetylene 
(C2H2)  as  a  soot  precursor  [140],  and  soot  oxidation  employing 
the  Nagle-Strickland-Constable  model  [141].  For  NOx  emissions, 
a  mechanism  derived  from  the  Gas  Research  Institute  (GRI)  NO 
mechanism  [84]  was  used. 

During  the  simulations,  the  charge  was  assumed  homogeneous 
at  intake  valve  closure.  The  researcher  generated  an  equivalence 
ratio  versus  peak  temperature  operating  map  based  on  numeri¬ 
cal  simulations,  and  compared  the  DME  operation  to  n-heptane 
operation.  High  CO  concentration  was  found  at  low  peak  combus¬ 
tion  temperatures  and  rich  mixtures.  CO  emissions  for  DME  were 
greater  than  those  for  n-heptane,  because  of  a  lower  initial  mole 
fraction  of  oxygen.  HC  emissions  for  DME  and  n-heptane  were  quite 
similar  at  lean  mixtures.  However,  HC  emissions  for  DME  were 
lower  than  those  for  n-heptane  below  a  peak  cycle  temperature 
of  1600K  and  an  equivalence  ratio  higher  than  1.5.  Although  such 
operating  regimes  are  outside  the  usual  HCCI  operating  conditions, 
this  finding  indicated  that  at  low  temperature-high  equivalence 
ratio  combustion  conditions  HC  emissions  can  be  reduced  by  sub¬ 
stituting  DME  for  diesel  fuel.  The  optimal  operating  range  for  the 
best  ISFC  or  imep  is  at  low  peak  cycle  temperature  and  below  stoi¬ 
chiometric  for  both  fuels,  because  ignition  is  near  TDC  for  that  area. 
Total  heat  release  of  DME  is  close  to  that  of  n-heptane  in  this  study. 
The  optimal  operating  range,  which  was  based  on  a  merit  func¬ 
tion  including  emissions  and  performance,  was  between  1700  and 
1850  K  peak  cycle  temperature  and  below  0.5  equivalence  ratio  for 
DME  HCCI  engines. 

Kong  [142]  studied  numerically  and  experimentally  the  dual 
fuel  DME-methane  oxidation  in  an  HCCI  engine.  The  researcher 
used  an  improved  version  of  the  CFD  code  KIVA-3V  [121].  The 
CHEMKIN  chemistry  solver  [94]  was  integrated  into  the  KIVA-3V 
code  for  solving  the  chemistry  during  multi-dimensional  engine 
simulations.  A  detailed  reaction  mechanism  for  DME  [53-55]  was 
used  to  simulate  fuel  chemistry,  consisting  of  79  species  and  351 
reactions.  This  mechanism  was  enriched  by  NOx-related  reactions 
obtained  from  a  NO  mechanism,  which  was  derived  from  the  Gas 
Research  Institute  (GRI)  NO  mechanism  [84].  The  final  scheme  con¬ 
sisted  of  83  species  and  360  reactions.  These  reactions  included  the 
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oxidation  chemistry  of  methane,  which  was  used  to  simulate  nat¬ 
ural  gas.  Turbulence  was  modeled  via  the  RNG  k-s  model  and  a 
new  wall  function  was  constructed  to  calculate  the  wall  heat  flux, 
which  accounted  for  the  variation  of  the  gas  density  and  the  turbu¬ 
lent  Prandtl  number  in  the  boundary  layer  assuming  quasi-steady 
conditions. 

Using  this  model,  the  researcher  found  that  as  DME  concen¬ 
tration  decreased  ignition  delayed  and  combustion  duration  was 
prolonged.  In  contrast  to  methane,  DME  experienced  two-stage 
ignition  and  was  depleted  before  any  significant  oxidation  of 
methane.  During  the  first  stage  of  heat  release,  a  noticeable  amount 
of  CO  was  formed,  which  was  later  increased  during  combustion 
and  decreased  rapidly  only  after  both  methane  and  DME  were 
nearly  consumed.  The  combustion  of  methane  was  induced  by  the 
high  temperatures  resulting  from  DME  combustion.  NO*  emissions 
were  much  lower  than  those  in  conventional  diesel  engines,  due 
to  the  low  combustion  temperature.  It  was  found  that  there  was  a 
minimum  requirement  in  DME  concentrations  for  HCCI  combus¬ 
tion  regardless  of  natural  gas  concentrations,  in  order  to  obtain 
stable  combustion.  Moreover,  the  operating  limits  -  defined  by 
misfire  and  knocking  combustion  -  narrowed  as  the  methane  con¬ 
centration  increased,  implying  that  high  methane  concentrations 
induce  combustion  instability  under  the  conditions  examined.  This 
was  probably  the  result  of  the  rapid  combustion  of  methane,  which 
was  triggered  by  the  DME  oxidation  and  prevailed  as  the  amount 
of  methane  increased. 

Yu  et  al.  [143]  used  large  eddy  simulation  (LES)  to  simulate 
ethanol  HCCI  combustion,  in  which  the  ethanol  was  port  injected. 
The  main  scope  of  the  study  was  to  investigate  the  role  of  turbu¬ 
lence  in  the  HCCI  combustion  process.  The  chemical  mechanism 
used  was  that  of  Marinov  [60].  The  LES  simulations  were  conducted 
for  two  different  piston  configurations,  i.e.  a  disc  shaped  combus¬ 
tion  chamber  and  a  bowl  in  piston  configuration,  while  maintaining 
the  compression  ratio.  A  grid  consisting  of  128  x  128  x  128  mesh 
cells  was  used,  with  size  about  1  mm  in  the  radial  direction  and 
0.32-1.5  mm  in  the  axial  direction  for  the  square  bowl-in-piston 
case  and  0.075-1.3  mm  for  the  disc-shaped  combustion  chamber 
case.  This  mesh  size  was  of  the  order  of  Taylor  microscale  and  at 
TDC  close  to  the  Kolmogorov  microscale.  Eight  parallel  processors 
and  48  h  were  necessary  to  yield  a  full  cycle  simulation  for  a  four- 
stroke  engine.  The  simulation  results,  which  were  validated  against 
experiments,  showed  that  the  piston  geometry  affected  the  HCCI 
combustion  process;  the  square  bowl-in-piston  design  had  a  much 
slower  combustion.  Although  the  levels  of  turbulence  in  terms  of 
the  velocity  components  were  similar,  the  temperature  inhomo¬ 
geneities  were  quite  different  between  the  two  configurations.  This 
was  the  reason  for  the  longer  combustion  duration  in  the  case  of 
square  bowl-in-piston  design. 

A  similar  study  was  later  conducted  by  Joelsson  et  al.  [144].  In 
this  study,  the  researchers  compared  numerically  and  experimen¬ 
tally  the  ethanol  HCCI  combustion  in  a  metal  and  an  optical  engine. 
The  engines  were  essentially  identical,  the  main  difference  being 
the  piston  material;  a  metal  piston  was  used  in  the  metal  engine 
and  a  quartz  piston  head  in  the  optical  engine.  It  was  expected  that 
the  quartz  piston  engine  would  operate  at  a  higher  piston  tempera¬ 
ture,  since  the  quartz  piston  was  air  cooled,  in  contrast  to  the  metal 
piston  that  was  oil  cooled.  The  scope  of  the  study  was  to  investi¬ 
gate  the  resulting  thermal  stratification  differences  and  their  effect 
on  ethanol  HCCI  combustion.  The  same  LES  model  was  used  in  this 
study,  as  in  the  previous  one.  The  expected  difference  in  the  piston 
temperatures  was  accounted  for  by  increasing  the  piston  surface 
temperature  in  the  optical  HCCI  engine  simulation.  The  simulation 
results  showed  that  the  two  engines  had  minor  differences  in  the 
induced  turbulence  field.  However,  the  thermal  stratification  was 
much  higher  in  the  optical  engine,  due  to  the  hotter  quartz  piston 
temperature  and  the  lower  intake  temperature  used  in  the  optical 


engine.  These  factors  prolonged  the  combustion  duration,  thereby 
decreasing  the  maximum  pressure  rise  rate  and  the  peak  pressure. 

Viggiano  and  Magi  [145]  used  a  multidimensional  CFD  model 
coupled  with  chemical  kinetics  to  study  ethanol  HCCI  combustion. 
The  model  used  was  the  REC-200,  which  solves  the  3D  Reynolds 
averaged  Navier-Stokes  equations  for  transient,  two-phase,  tur¬ 
bulent,  chemically  reactive  flows  with  sprays  [146].  Turbulence 
was  modeled  with  a  two-equation  k-s  model,  and  the  heat  and 
momentum  fluxes  were  calculated  using  the  wall  function  model  of 
Launder  and  Spalding  [147].  Combustion  was  modeled  by  coupling 
the  turbulent  flow  to  a  detailed  reaction  mechanism.  The  approach 
used  for  the  coupling  was  that  of  Kong  and  Reitz  [148].  The  influence 
of  turbulence  on  the  reaction  rates  was  accounted  for  by  multiply¬ 
ing  the  latter  with  a  correction  factor,  which  accounted  for  both 
the  kinetic  and  turbulence  timescales.  The  detailed  reaction  mech¬ 
anism  used  was  that  of  Saxena  and  Williams  [59],  which  included 
43  species  and  235  reactions  excluding  NOx  formation  reactions. 
The  structured  grid  of  the  1 -degree  sector  considered,  consisted 
of  26  cells  in  the  radial  direction.  In  the  axial  direction,  the  cell 
number  varied  from  48  at  BDC  to  8  at  TDC  by  removing  cell  rows 
during  compression.  The  grid  resolution  near  the  wall  was  kept 
at  0.2  mm.  After  the  validation  of  the  CFD  model  against  experi¬ 
ments  was  accomplished,  the  researchers  investigated  the  effect 
of  turbulent  timescale,  turbulent  diffusivity,  initial  and  wall  tem¬ 
perature,  and  equivalence  ratio  on  HCCI  combustion.  They  found 
that  increasing  the  turbulent  timescale  induced  a  reduction  to  the 
combustion  rate  and  decreased  the  pressure  rise  rate.  Increasing 
the  initial  diffusivity  delayed  ignition  and  deteriorated  the  com¬ 
bustion  rate.  Inlet  and  wall  temperature  affected  combustion  in  a 
similar  manner,  with  the  effect  of  inlet  temperature  being  more 
pronounced;  specifically,  increasing  Tin  or  Tw  advanced  ignition, 
enhanced  combustion,  and  reduced  CO  emissions. 

In  Ref.  [149]  Liu  and  Karim  examined  numerically  the  effect  of 
various  operating  parameters  on  hydrogen  HCCI  combustion.  For 
the  purposes  of  this  investigation,  the  KIVA3  multi-dimension  CFD 
code  [121]  was  adopted  with  some  modifications  to  make  it  suit¬ 
able  for  modeling  the  turbulent  combustion  processes.  Detailed 
chemical  kinetics  of  hydrogen  oxidation  was  implemented  into 
the  code,  including  20  reactions  with  additional  reactions  involv¬ 
ing  nitrogen  [71,72].  With  this  model,  the  researchers  conducted 
a  parametric  study,  the  results  of  which  showed  that  by  increas¬ 
ing  the  swirl  ratio  heat  transfer  was  enhanced,  thereby  reducing 
both  the  minimum  and  maximum  temperature  within  the  com¬ 
bustion  chamber.  As  a  result,  ignition  was  retarded,  combustion 
duration  was  prolonged,  and  the  peak  pressure  and  NOx  emissions 
were  reduced.  An  increase  in  engine  speed  increased  the  heat  loss 
rate  due  to  the  induced  increase  of  the  swirl  ratio,  but  the  heat 
losses  were  reduced  because  the  time  available  for  heat  transfer 
was  decreased.  The  calculations  showed  that  there  was  an  opti¬ 
mum  engine  speed  in  terms  of  indicated  thermal  efficiency.  Initial 
temperature  and  compression  ratio  affected  significantly  the  onset 
of  combustion  as  well  as  the  combustion  duration. 

Cotino  and  Jeanmart  [150]  studied  numerically  the  running  zone 
of  ethyl  acetate  and  compare  it  with  the  running  zone  of  isooctane. 
Ethyl  acetate  was  chosen,  since  it  is  a  product  of  the  acidogenic 
fermentation  process  used  to  produce  biodiesel.  This  molecule  was 
the  basic  compound  and  was  used  as  a  starting  point  to  character¬ 
ize  a  group  of  esters  produced.  The  model  used  in  the  numerical 
study  was  the  commercial  CFD  FLUENT  code  and  the  computa¬ 
tional  grid  was  obtained  from  the  commercial  software  GAMBIT. 
The  combustion  mechanism  used  for  ethyl  acetate  was  a  reduced 
one,  consisting  of  41  species  and  203  reactions.  This  mechanism 
was  derived  from  more  detailed  ones,  namely  the  Gasnot  et  al. 
ethyl  acetate  submechanism  [65],  and  the  low  carbon  combustion 
mechanism  of  Miller  and  Melius  [151].  Turbulence  was  modeled 
with  the  RNG  k-s  model  and  the  molecular  mixing  at  the  subgrid 
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level  via  the  eddy  dissipation  model  (EDC).  Except  for  swirl  motion, 
3D  effects  were  not  included  and  an  axisymmetric  geometry  was 
chosen,  yielding  33,000  cells  at  IVC. 

The  model  was  validated  against  experimental  results  using 
isooctane  as  the  fuel,  and  was  subsequently  used  with  ethyl  acetate 
to  determine  the  ethyl  acetate  HCCI  operating  limits  (running  zone) 
at  1 200  rpm.  The  upper  load  limit  was  determined  by  the  maximum 
pressure  rise  rate,  which  was  set  at  5  bar/CAD.  The  lower  limit  cri¬ 
terion  implied  combustion  instability  and  was  determined  taking 
into  account  the  combustion  timing,  the  EGR  rate,  the  sensitivity  of 
the  combustion  on  the  initial  temperature,  the  compatibility  of  the 
final  and  initial  temperatures,  and  the  CO  emissions.  The  results  of 
the  study  showed  that  the  ethyl  acetate  had  a  similar  running  zone 
range  as  isooctane,  but  it  was  shifted  toward  higher  loads.  More¬ 
over,  ethyl  acetate  presented  extended  low  temperature  reactions 
compared  to  isooctane. 

Um  and  Park  [152]  analyzed  the  chemical  effects  of  biodiesel 
on  combustion  and  emissions  of  an  HCCI  engine  fueled  with 
diesel/biodiesel  blends.  For  this  purpose,  they  used  a  modified  ver¬ 
sion  of  the  KIVA-3  V  code  [  1 2 1  ]  in  a  2D  computational  mesh,  coupled 
with  CHEMKIN  chemistry  solver  [70].  The  mechanism  for  biodiesel 
oxidation  used  was  the  one  suggested  by  Brakora  et  al.  [78],  and 
was  modified  by  Um  and  Park  [153]  to  produce  a  reasonable  igni¬ 
tion  delay.  The  mechanism  consisted  of  156  reactions  among  54 
species.  Since  it  was  found  that  neat  methyl  butanoate  was  not 
suitable  as  a  biodiesel  surrogate,  due  to  its  high  oxygen  content,  it 
was  assumed  that  1  mol  of  biodiesel  consisted  of  1  mol  of  methyl 
butanoate  and  2  mol  of  n-heptane.  Neat  heptane  was  also  used  as  a 
diesel  surrogate.  A  NOx  mechanism,  reduced  from  the  Gas  Research 
Institute  (GRI)  NO  mechanism  [84],  was  used  for  this  study.  It  con¬ 
sisted  of  12  reactions  and  4  species  (i.e.,  N,  NO,  N02,  and  N20). 
Moreover,  a  soot  formation  and  oxidation  model  was  incorporated, 
since  rich  mixtures  with  equivalence  ratio  up  to  2.0  were  used.  The 
high  equivalence  ratio,  which  is  impractical  in  actual  HCCI  engines, 
is  justified  considering  the  nature  of  the  study,  which  was  to  exam¬ 
ine  the  chemical  properties  of  the  biodiesel  blends.  Soot  formation 
was  calculated  using  a  modified  Hiroyasu  model,  which  considers 
acetylene  (C2H2)  as  a  soot  precursor  [140],  while  soot  oxidation 
was  modeled  using  the  Nagle-Strickland-Constable  model  [141]. 

The  results  of  the  study  revealed  that  at  stoichiometric  equiva¬ 
lence  ratio,  increasing  the  biodiesel-to-diesel  fuel  ratio  advanced 
ignition  and  decreased  the  total  amount  of  heat  released.  This 
reduced  the  in-cylinder  temperatures  and  induced  a  decrease  in 
NO*  emissions  with  an  increase  in  CO  and  HC.  At  richer  mixtures 
(0>1.O)  increasing  the  biodiesel  content  increased  CO,  decreased 
NO*  and  HC  emissions,  and  decreased  Soot  emissions  compared  to 
neat  diesel  fuel  operation.  Below  stoichiometric  conditions  (0  <  1 .0) 
CO  and  HC  were  formed  in  the  cases  where  peak  cycle  temperatures 
were  below  1 400  K,  NOx  was  produced  at  peak  temperatures  higher 
than  1800K,  while  Soot  emissions  were  minimal.  At  rich  mixtures 
significant  amounts  of  CO  and  HC  emissions  were  observed  for  most 
cases  and  soot  emissions,  which  were  higher  for  neat  diesel,  peak¬ 
ing  in  the  cases  with  a  peak  cycle  temperature  of  1600K.  The  NOx 
and  CO  emission  maps,  i.e.  equivalence  ratio  versus  peak  tempera¬ 
ture  maps,  were  similar  for  the  diesel  fuel  and  biodiesel  cases. 


6.  Conclusions 

Chemical  kinetics  are  an  essential  part  of  biofueled  HCCI  com¬ 
bustion  modeling,  owing  to  the  nature  of  HCCI  combustion  and 
the  diversity  of  the  biofuels’  molecular  structure.  The  development 
of  chemical  reaction  mechanisms  is  an  active  research  area,  from 
which  biofuel  HCCI  research  can  benefit.  Efforts  should  be  focused 
on  the  consistency  of  the  various  mechanisms  developed  for  use  in 
HCCI  engines.  Special  care  should  be  exercised  on  the  selection  of 


a  chemical  mechanism,  depending  on  the  conditions  under  which 
its  validation  was  conducted.  These  conditions  are  preferably 
high  pressure  and  low  equivalence  ratio,  since  such  conditions 
are  encountered  in  HCCI  combustion.  It  would  also  be  advisable 
to  implement  more  than  one  chemical  mechanism  for  the  same 
biofuel  in  any  HCCI  simulation  model,  and  compare  the  results. 
This  helps  to  avoid  reaching  conclusions  (e.g.  on  autoignition,  or 
emissions  formation),  which  are  the  result  of  the  specific  chemical 
mechanism  used  rather  than  the  HCCI  combustion  process. 

Single-zone  models  for  biofuel  HCCI  combustion  are  very  use¬ 
ful  in  identifying  reaction  paths  and  dominant  chemical  species 
during  combustion.  They  have  been  used  to  provide  insight  into 
the  behavior  of  various  biofuels  and  blends  thereof,  and  to  validate 
the  corresponding  combustion  mechanisms  versus  experiments  or 
other  more  detailed  mechanisms.  Owing  to  the  zero  dimensional 
nature  of  the  single-zone  models,  they  have  been  used  mainly  to 
provide  insight  into  the  autoignition  properties  of  the  fuel  and  the 
formation  of  NOx  emissions.  Single-zone  models  also  provide  an 
indication  on  the  effect  of  various  configurations,  such  as  nega¬ 
tive  valve  overlap,  heavy  EGR,  fuel  reforming,  etc.,  on  biofuel  HCCI 
combustion. 

Multi-zone  models  with  their  added  feature  of  accounting  for 
mixture  thermal  and  fuel  stratification,  provide  better  and  more 
realistic  estimates  on  combustion  duration,  peak  combustion  pres¬ 
sure  and  emissions  formation.  Multi-zone  models  for  biofuel  HCCI 
combustion  have  been  used  to  study  the  effect  of  various  biofuels 
and  their  blends  on  the  in-cylinder  thermal  stratification,  combus¬ 
tion  duration,  emissions  formation,  and  the  overall  performance. 
Not  only  can  they  provide  an  insight  into  the  formation  of  CO  and 
HC  emissions,  but  also  to  the  specific  species  formed  and  consumed 
in  the  low  temperature  regions  of  the  combustion  chamber,  thereby 
indicating  the  species  expected  to  be  found  in  the  exhaust  of  a  HCCI 
engine.  This  is  an  important  feature,  since  the  composition  of  HC 
emissions  can  vary  substantially  in  the  case  of  biofuels,  due  to  their 
non-conventional  molecular  structure. 

Multi-dimensional  models  provide  an  even  higher  spatial  res¬ 
olution  of  the  combustion  chamber  than  the  multi-zone  ones.  Gas 
flow  during  the  gas  exchange  process  or  within  the  combustion 
chamber  can  be  modeled  in  more  detail,  providing  a  more  realis¬ 
tic  spatial  resolution  of  the  species  evolution  during  combustion. 
Moreover,  mass  flow  between  hot  and  cold  regions  is  described  in 
detail,  providing  an  insight  into  the  emissions  formation  process. 
Especially  in  the  case  of  direct  injection  biofueled  HCCI  combus¬ 
tion,  the  physical  properties  of  the  fuel  must  be  described  in  detail; 
properties  such  as  viscosity,  volatility,  etc.,  are  essential  to  describe 
the  evolution  of  the  fuel  spray  into  the  combustion  chamber.  The 
importance  of  these  properties  cannot  be  overemphasized,  since 
the  spray  evolution  will  determine  not  only  the  fuel  stratification 
within  the  combustion  chamber  but  also  the  fuel  impingement  on 
the  cylinder  liner  or  the  piston  crown.  The  description  of  the  physic¬ 
ochemical  properties  of  biofuels  is  an  active  area  of  research,  which 
is  bound  to  be  further  advanced. 
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